Fish and health with national fish promotion programme
This is an open assessment looking at potential health effects of a national fish promotion program in Finland. The details of the assessment are described at Opasnet. This file contains the R code to run the assessment model.
Knit to html for best performance.
Calculation is based on BONUS GOHERR project and its http://en.opasnet.org/w/Goherr_assessment.
What needs to be done for PFAS assessment?
- Amount should be gender and age-specific, because we need to worry about young mothers. Solution: take KKE amounts of fish, and scale those with the Goherr subgroup-specific proportions. DONE
- Amounts should reflect the actual fish consumption in Porvoo. Solution: postpone and use national statistics. DONE
- Infant’s dioxin concentration module must be added. Solution: Use Goherr model for expo_indir DONE
- The module must be updated to match PFAS as well. Solution: update the module to contain column Exposure_agent. Change body fat parameter to volume of distribution.
- Find why expo_indir is so much higher than EFSA toxicokinetic assessment. Make an alternative model.
- Add ERF for PFAS (sum of PFOS, PFHxS, PFOA, PFNA). Solution: Make a new page for ERF of PFAS and add that to adjusted ERFs. First immunology; postpone cholesterol and low birth weight and liver toxicity. DONE
- Add case burdens for PFAS outcomes. Solution: give rough estimates for immunology to get started. DONE
- Add PFAS concentrations to data. Solution: Look at Porvoo concentrations first and make conc_pfas; combine that with conc_vit. DONE
- Add dioxin and MeHg concentrations. Later.
- Now that the model runs technically, look through each part to check that it makes sense.
- If PAF contains some inputs with Age and some without Age (implying that they apply to any age group), only the explicit age groups remain in the end. Add RRorig to compensate for this. DONE
- TEQ and MeHg are missing from dose. Check why and correct.
- The proportions of All-cause mortality and CHD2 mortality are 3:1 in http://fi.opasnet.org/fi/Kotimaisen_kalan_edist%C3%A4misohjelma but 3:2 in this analysis. Why?
- In exposure, ALA + DHA < Omega3. Why?
- BoDattr: age group 45 - 50 is way smaller than any other group. Why?
- conc_mehg does not contain Iter although it is probabilistic. Update formula.
objects.latest("Op_en2261", code_name="RR2") # RR on page [[Health impact assessment]]## Loading objects:
## RRRRorig <- RR@formula
RR@formula <- function(...) {
out <- RRorig()
out <- out * Ovariable(output=data.frame(Age=c(
"0 - 4", "5 - 9", "10 - 14", "15 - 19", "20 - 24", "25 - 29", "30 - 34", "35 - 39", "40 - 44", "45 - 49", "50 - 54",
"55 - 59", "60 - 64", "65 - 69", "70 - 74", "75 - 79", "80 - 84", "85+", "Undefined"),
Result=1),
marginal=c(TRUE,FALSE)
)
return(out)
}# This code was forked from https://github.com/jtuomist/fishhealth/blob/master/fishhealth.Rmd
# This code was previously forked from code Op_fi5923/model on page [[Kotimaisen kalan edistämisohjelma]]
# The code was even more previously forked from Op_fi5889/model on page [[Ruori]] and Op_en7748/model on page [[Goherr assessment]]
dat <- opbase.data("Op_fi5932", subset="Malliparametrit")[-1] # [[PFAS-yhdeisteiden tautitaakka]]
dec <- opbase.data("Op_fi5932", subset="Decisions")[-1]
DecisionTableParser(dec)
CTable <- opbase.data("Op_fi5932",subset="CollapseMarginals")
#for(i in 1:ncol(CTable)) {CTable[[i]] <- as.character(CTable[[i]])} # The default is currently character, not factor
CollapseTableParser(CTable)
cat("Laskennassa käytetty data.\n")## Laskennassa käytetty data.datcat("Tarkastellut päätökset.\n")## Tarkastellut päätökset.deccat("Aggregoidut marginaalit.\n")## Aggregoidut marginaalit.CTabledummy <- Ovariable("dummy",data=data.frame(Age="dummy", Fish="dummy", Compound="dummy", Area="dummy", Result=1)) # Keep these columns marginals
fish_proportion <- Ovariable( # How population subgroups eat fish differently
"fish_proportion",
data = data.frame(prepare(dat,"fish proportion",c("Type","Exposure_agent","Response","Unit"))),
unit="proportion of the mean")
total_amount <- Ovariable(
"total_amount",
data=prepare(dat, "amount", c("Type","Response","Exposure_agent","Unit")),
unit="M kg/a")
amount <- Ovariable(
"amount",
dependencies = data.frame(Name=c("total_amount", "fish_proportion")),
formula = function(...) {
amount <- total_amount
# Filleted weight, i.e. no loss.
amount <- amount * 1000 / 5.52 /365.25
# M kg/a per 5.52M population --> g/d per average person.
amount <- amount * fish_proportion
# fish_proportion tells the relative amount in each subgroup
# Match KKE-classification in amount with Fineli classification
tmp <- Ovariable(
output = data.frame(
Kala = c("Kasvatettu", "Kaupallinen", "Kirjolohi", "Silakka", "Vapaa-ajan", "Muu tuonti", "Tuontikirjolohi", "Tuontilohi"),
Fish = c("Whitefish", "Average fish","Rainbow trout", "Herring", "Average fish", "Average fish", "Rainbow trout", "Salmon"),
Result = 1
),
marginal = c(TRUE, TRUE, FALSE)
)
amount <- amount * tmp
return(amount)
},
unit="g/d"
)
# Exposure:To child and To eater not needed, because dioxins are not (yet) included
population <- Ovariable(
"population",
data=prepare(dat, "population", c("Type", "Exposure_agent", "Response","Unit")),
unit = "#"
)
population@data$Age <- factor(population@data$Age, levels=c(
"0 - 4", "5 - 9", "10 - 14", "15 - 19", "20 - 24", "25 - 29", "30 - 34", "35 - 39", "40 - 44", "45 - 49", "50 - 54",
"55 - 59", "60 - 64", "65 - 69", "70 - 74", "75 - 79", "80 - 84", "85+", "Undefined"))
incidence <- Ovariable(
"incidence",
data = prepare(dat,"incidence",c("Type","Exposure_agent","Unit")),
unit="1/person-year")
#incidence@data$Age[is.na(incidence@data$Age)] <- ""
case_burden <- Ovariable(
"case_burden",
data = prepare(dat,"case burden",c("Type", "Exposure_agent","Unit")),
unit="DALY/case")
ERFchoice <- Ovariable(
"ERFchoice",
data =
prepare(dat, "ERFchoice", c("Unit", "Type"))
)
InpBoD <- EvalOutput(Ovariable( # Evaluated because is not a dependency but an Input
"InpBoD",
data = prepare(dat, "BoD", c("Type","Exposure_agent","Unit")),
unit="DALY/a"
))
InpBoD$Response[InpBoD$Response=="All causes"] <- "All-cause mortality"
InpBoD$Response[InpBoD$Response=="Depressive disorders"] <- "Depression"
InpBoD$Response[InpBoD$Response=="Neoplasms"] <- "Cancer morbidity"
InpBoD$Response[InpBoD$Response=="Respiratory infections and tuberculosis"] <- "Immunosuppression" # Infections of 0-9-year-olds are assumed to represent the background BoD of immunosuppressive effect of PFAS
InpBoD$Response[InpBoD$Response=="Cardiovascular diseases"] <- "CHD2 mortality"conc_vit <- Ovariable(
"conc_vit",
ddata = "Op_en1838", # [[Concentrations of beneficial nutrients in fish]]
subset = "Fineli data for common fish species",
unit="ALA mg/g, DHA mg/g, Fish g/g, Omega3 mg/g, Vitamin D µg/g f.w. after adjustment"
)
df = conc_vit@data
df$Nutrient[df$Nutrient=="D-vitamiini (µg)"] <- "Vitamin D"
df$Nutrient[df$Nutrient=="rasvahapot n-3 moni-tyydyttymättömät (g)"] <- "Omega3"
df$Nutrient[df$Nutrient=="rasvahappo 18:3 n-3 (alfalinoleenihappo) (mg)"] <- "ALA"
df$Nutrient[df$Nutrient=="rasvahappo 22:6 n-3 (DHA) (mg)"] <- "DHA"
df$Nutrient[df$Nutrient=="proteiini (g)"] <- "Fish"
df$conc_vitResult[df$Nutrient=="Fish"] <- "1"
df <- dropall(df[df$Nutrient %in% c("Vitamin D", "Omega3", "ALA", "DHA", "Fish") , ])
conc_vit@data <- df
######## Concentration of PFAS
# Data from EU-kalat3 (Finland excl Vanhankaupunginlahti): # pg/g fresh weight
# POP mean sd min Q0.025 median Q0.975 max
# 2.5% PFOS 2055.757 1404.045 305.0399 330.1365 1533.269 5029.697 5814.935
# Data from EU-kalat3 (Vanhankaupunginlahti, Helsinki) # ng/g f.w.
# POP mean sd min Q0.025 median Q0.975 max
#2.5% PFOS 14.428 11.94542 1.499441 1.607789 15.64988 35.32517 38.91994
conc_eukalat <- EvalOutput(Ovariable(
"conc_eukalat",
data = data.frame(
Area = c("Suomi","Helsinki"),
Compound="PFOS",
Result=c("2.056 (3.301 - 5.030)", "14.428 (1.499 - 35.325)")),
unit="ng/g fresh weight"
))
conc_pfas_raw <- Ovariable(
"conc_pfas_raw",
ddata="Op_fi5932",subset="PFAS concentrations",
unit="ng/g f.w.")
conc_pfas <- Ovariable(
"conc_pfas",
dependencies = data.frame(Name="conc_pfas_raw"),
formula = function(...) {
out <- conc_pfas_raw
out$Source <- NULL
out$Area <- "Porvoo"
out <- expand_index(out, list(Fish=list(Perch=c(
"Average fish", "Pike","Rainbow trout","Roach", "Salmon", "Vendace", "Whitefish"))))
out@marginal <- c(TRUE, TRUE, FALSE, TRUE, TRUE)
return(out)
},
unit="ng/g fresh weight"
)
#sum_pfas <- oapply(conc_pfas, cols=c("Kala","Compound"), FUN=sum)
#tmp <- conc_pfas / sum_pfas
#summary(tmp, marginals="Compound")
#
## This tells that PFOS consists of 71 - 97 % of the four key PFAS, while PFOA, PFNA, and PFHxS consist of
# 0 - 10 %, 2 - 18 %, and 0 - 9 %, respectively.
# Even if we included the next most abundant congeners, i.e. PFDA and PFUnA, the overall picture would not change.
conc <- Ovariable(
"conc",
dependencies = data.frame(
Name=c("conc_vit", "conc_pfas","conc_pcddf","conc_mehg"),
Ident = c(NA,NA,"Op_en3104/pcddf_mean","Op_en4004/conc_mehg_allfish")),
formula = function(...){
conc_vit <- oapply(conc_vit, cols=c("Kala", "Adjust"),FUN=mean)
colnames(conc_vit@output)[colnames(conc_vit@output)=="Nutrient"] <- "Compound"
conc_pfas <- oapply(conc_pfas, cols=c("Obs","Area"), FUN=mean)
conc_pfas$Compound[conc_pfas$Compound %in% c("PFOA","PFNA","PFHxS","PFOS")] <- "PFAS"
conc_pfas <- oapply(conc_pfas, cols="", FUN=sum)
conc_pcddf$Compound <- "TEQ"
conc_mehg <- oapply(conc_mehg[conc_mehg$Area=="Porvoo",], NULL, mean, c("Area","Kala")) # UPDATE FOR PROBABILISTIC APPROACH
conc_mehg$Compound <- "MeHg"
out <- OpasnetUtils::combine(conc_vit, conc_pfas, conc_pcddf, conc_mehg)
return(out)
}
)
#conc <- EvalOutput(conc, verbose=TRUE)
#View(conc@output)
#View(conc_pcddf@output)
#View(conc_mehg@output)###################################################################
# Code copied from http://en.opasnet.org/w/Goherr_assessment#
mc2dparam<- list(
N2 = 10, # Number of iterations in the new Iter
strength = 50, # Sample size to which the fun is to be applied. Resembles number of observations
run2d = FALSE, # Should the mc2d function be used or not?
info = 1, # Ovariable that contains additional indices, e.g. newmarginals.
newmarginals = c("Exposure"), # Names of columns that are non-marginals but should be sampled enough to become marginals
method = "bootstrap", # which method to use for 2D Monte Carlo? Currently bootsrap is the only option.
fun = mean # Function for aggregating the first Iter dimension.
)
if(FALSE) {
## Exposure with background exposure but without mother's exposure to child
expo_dir <- Ovariable(
"expo_dir",
dependencies=data.frame(Name=c("amount","conc","expo_bg")),
formula = function(...) {
out <- conc[conc$Exposure_agent=="TEQ",] * 0 + 1
out$Exposure_agent <- "Fish"
out <- combine(conc, out, name="conc")
out <- oapply(amount * out, cols="Fish", FUN=sum)
out <- Ovariable(output = data.frame(
Exposcen = c("BAU", "No exposure"),
Result = c(1, 0)
), marginal=c(TRUE,FALSE)) * out + expo_bg
out$Exposure <- as.factor(
ifelse(
out$Exposure_agent %in% c("DHA", "MeHg"),
"To child",
"To eater"
)
)
return(out)
},
unit = "PCDDF, PCB, TEQ: pg /d; Vitamin D, MeHg: µg /d; DHA, EPA, Omega3: mg /d; Fish: g /d"
)
## Background-exposure to vitamin D and omega-3
addexposure <- Ovariable(
"addexposure",
ddata = "Op_en7748", # [[Benefit-risk assessment of Baltic herring and salmon intake]]
subset = "Background exposure",
unit = "PCDDF, PCB, TEQ: pg /d; Vitamin D, MeHg: µg /d; DHA, EPA, Omega3: mg /d"
)
# Should the background be specific for gender and country? At the moment it is.
expo_bg <- Ovariable(
"expo_bg",
dependencies = data.frame(Name="addexposure","info"),
formula = function(...) {
out <- addexposure
# Empty values ("") in indices must be replaced by NA so that Ops works correctly.
levels(out$Gender)[levels(out$Gender) == ""] <- NA
levels(out$Country)[levels(out$Country) == ""] <- NA
levels(out$Exposure_agent)[levels(out$Exposure_agent) == ""] <- NA
out@output <- fillna(out@output, c("Country", "Gender", "Exposure_agent"))
temp1 <- out[out$Exposure_agent %in% c("PCDDF","PCB") , ]
temp1 <- oapply(temp1, cols = "Exposure_agent", FUN = sum)
temp1$Exposure_agent <- "TEQ"
temp2 <- out[out$Exposure_agent %in% c("EPA", "DHA") , ]
temp2 <- oapply(temp2, cols = "Exposure_agent", FUN = sum)
temp2$Exposure_agent <- "Omega3"
out <- OpasnetUtils::combine(out, temp1, temp2)
out <- unkeep(out * info, prevresults = TRUE, sources = TRUE)
return(out)
},
unit = "PCDDF, PCB, TEQ: pg /d; Vitamin D, MeHg: µg /d; DHA, EPA, Omega3: mg /d"
)
# Stores non-marginal columns for further use.
info <- Ovariable(
"info",
dependencies = data.frame(Name = c("jsp")),
formula = function(...) {
out <- jsp
out$Group <- factor(
paste(out$Gender, out$Ages),
levels = c("Female 18-45", "Male 18-45", "Female >45", "Male >45")
)
out$Country <- factor(out$Country, ordered=FALSE)
out <- unique(out@output[c("Iter","Country","Group","Gender","Row")])
out$Result <- 1
return(out)
}
)
} # END IF
############################### Code from Goherr assessment ends
expo_dir <- Ovariable(
"expo_dir",
dependencies = data.frame(Name=c("conc", "amount")),
formula = function(...) {
conc$Fish[conc$Fish %in% c("Perch","Pike-perch","Eel","")] <- "Average fish"
conc <- oapply(conc, cols="", FUN=mean)
out <- conc * amount
# colnames(out@output)[colnames(out@output)=="Nutrient"] <- "Exposure_agent"
return(out)
},
unit="ng/d"
)
exposure <- Ovariable(
"exposure",
dependencies = data.frame(
Name = c(
"expo_dir", # direct exposure, i.e. the person eats or breaths the exposure agent themself
"expo_indir", # indirect exposure, i.e. the person (typically fetus or infant) is exposed via someone else (mother)
"mc2d" # 2D Monte Carlo function
),
Ident = c(
NA,
"Op_en7797/expo_indir2", # [[Infant's dioxin exposure]] # expo_indir
"Op_en7805/mc2d") # [[Two-dimensional Monte Carlo]]
),
formula = function(...) {
out <- OpasnetUtils::combine(expo_dir, expo_indir)
out <- unkeep(out, "Source.1", sources=TRUE)
out <- mc2d(out)
out$Exposure[is.na(out$Exposure)] <- "Direct"
return(out)
},
unit = "PCDDF, PCB, TEQ: (To eater: pg /day; to child: pg /g fat); Vitamin D, MeHg: µg /day; DHA, EPA, Omega3: mg /day"
)
exposure <- Ovariable(
"exposure",
dependencies=data.frame(Name=c("amount","conc")),
formula = function(...) {
out <- amount * conc
colnames(out@output)[colnames(out@output)=="Compound"] <- "Exposure_agent"
return(out)
}
)objects.latest("Op_en2261",code_name="BoDattr2") # [[Health impact assessment]]## Loading objects:
## BoDattrtryCatch(BoDattr <- EvalOutput(BoDattr, verbose=TRUE))## Evaluating BoDattr ...
##
## - BoD fetched successfully!
##
## - PAF fetched successfully!
## - Evaluating BoD ...
## - - Evaluating incidence ...
##
## done(0.01 secs)!
## - - Checking incidence marginals ... Response, Age, incidenceSource recognized as marginal(s).## Loading required package: reshape2##
## Attaching package: 'reshape2'## The following object is masked from 'package:tidyr':
##
## smiths## - - Processing incidence decisions ... done!
## - - Evaluating case_burden ...
##
## done(0 secs)!
## - - Checking case_burden marginals ... Response, case_burdenSource recognized as marginal(s).
## - - Processing case_burden marginal collapses ... done!
## - - Evaluating population ...
##
## done(0 secs)!
## - - Checking population marginals ... Gender, Age, populationSource recognized as marginal(s).
##
## - done(0.3 secs)!
## - Checking BoD marginals ... Response, Age, incidenceSource, Adjust, Gender, populationSource, BoDSource recognized as marginal(s).
## - Processing BoD inputs ... done!
## - Processing BoD marginal collapses ...## Warning in oapply(variable, FUN = fun[[i]], cols = cols[[i]], na.rm = TRUE):
## While oapplying BoD, found NAs in indices: Adjust, InpBoDSource. They were
## automatically filled using fillna, which may result in a multiplied population.
## Please check your ovariable before using oapply.## done!
## - Evaluating PAF ...
##
## - - dose fetched successfully!
##
## - - ERF fetched successfully!
##
## - - frexposed fetched successfully!
##
## - - P_illness fetched successfully!
##
## - - sumExposcen fetched successfully!
##
## - - mc2d fetched successfully!
## - - Evaluating dose ...
##
## - - - BW fetched successfully!
## - - - Evaluating exposure ...
## - - - - Evaluating amount ...
## - - - - - Evaluating total_amount ...
##
## done(0 secs)!
## - - - - - Checking total_amount marginals ... Kala, Scenario, total_amountSource recognized as marginal(s).
## - - - - - Evaluating fish_proportion ...
##
## done(0 secs)!
## - - - - - Checking fish_proportion marginals ... Gender, fish_proportionSource recognized as marginal(s).
##
## ---- done(0.12 secs)!
## - - - - Checking amount marginals ... Kala, Scenario, total_amountSource, Gender, fish_proportionSource, Fish, amountSource recognized as marginal(s).
## - - - - Evaluating conc ...
##
## - - - - - conc_pcddf fetched successfully!
##
## - - - - - conc_mehg fetched successfully!
## - - - - - Evaluating conc_vit ...
##
## done(0.02 secs)!
## - - - - - Checking conc_vit marginals ... Kala, Fish, Nutrient, conc_vitSource recognized as marginal(s).
## - - - - - Processing conc_vit decisions ... done!
## - - - - - Evaluating conc_pfas ...
## - - - - - - Evaluating conc_pfas_raw ...
##
## done(0.02 secs)!
## - - - - - - Checking conc_pfas_raw marginals ... Fish, Compound, conc_pfas_rawSource recognized as marginal(s).
##
## ----- done(0.09 secs)!
## - - - - - Checking conc_pfas marginals ... Fish, Compound, conc_pfas_rawSource, Area, conc_pfasSource recognized as marginal(s).
## - - - - - Evaluating conc_pcddf ...
##
## - - - - - - TEF fetched successfully!
## - - - - - - Evaluating TEF ...
##
## - - - - - - - TEFversion fetched successfully!
## - - - - - - - Evaluating TEFraw ...
##
## done(0.01 secs)!
## - - - - - - - Checking TEFraw marginals ... Group, Compound, TEFversion, TEFrawSource recognized as marginal(s).
##
## ------ done(27.63 secs)!
## - - - - - - Checking TEF marginals ... Compound, Group, TEFversion, TEFrawSource, TEFSource recognized as marginal(s).
##
## ----- done(1.03 mins)!
## - - - - - Checking conc_pcddf marginals ... Fish, conc_pcddfSource recognized as marginal(s).
## - - - - - Evaluating conc_mehg ...
##
## - - - - - - Hg fetched successfully!
##
## ----- done(23.03 secs)!
## - - - - - Checking conc_mehg marginals ... Area, Kala, Fish, conc_mehgSource recognized as marginal(s).
##
## ---- done(2.22 mins)!
## - - - - Checking conc marginals ... Fish, Compound, conc_pfas_rawSource, concSource recognized as marginal(s).
##
## --- done(2.22 mins)!
## - - - Checking exposure marginals ... Kala, Scenario, total_amountSource, Gender, fish_proportionSource, Fish, amountSource, Exposure_agent, conc_pfas_rawSource, concSource, exposureSource recognized as marginal(s).
## - - - Processing exposure marginal collapses ...## Warning in oapply(variable, FUN = fun[[i]], cols = cols[[i]], na.rm = TRUE):
## While oapplying exposure, found NAs in indices: conc_pfas_rawSource. They were
## automatically filled using fillna, which may result in a multiplied population.
## Please check your ovariable before using oapply.## done!
## - - - Evaluating BW ...
##
## done(0 secs)!
## - - - Checking BW marginals ... BWSource recognized as marginal(s).
##
## -- done(2.69 mins)!
## - - Checking dose marginals ... Scenario, total_amountSource, Gender, fish_proportionSource, amountSource, Exposure_agent, conc_pfas_rawSource, concSource, Scaling, exposureSource, BWSource, doseSource recognized as marginal(s).
## - - Processing dose marginal collapses ... done!
## - - Evaluating ERF ...
##
## - - - ERF_env fetched successfully!
##
## - - - ERF_omega3 fetched successfully!
##
## - - - ERF_mehg fetched successfully!
##
## - - - ERF_diox fetched successfully!
##
## - - - ERF_vit fetched successfully!
##
## - - - ERF_micr fetched successfully!
##
## - - - ERF_pfas fetched successfully!
## - - - Evaluating ERF_env ...
##
## done(0.02 secs)!
## - - - Checking ERF_env marginals ... Exposure_agent, Response, Subgroup, Exposure, ER_function, Scaling, Exposure_unit, Observation, ERF_envSource recognized as marginal(s).
## - - - Evaluating ERF_omega3 ...
##
## done(0.01 secs)!
## - - - Checking ERF_omega3 marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_omega3Source recognized as marginal(s).
## - - - Evaluating ERF_mehg ...
##
## done(0 secs)!
## - - - Checking ERF_mehg marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_mehgSource recognized as marginal(s).
## - - - Evaluating ERF_diox ...
##
## done(0 secs)!
## - - - Checking ERF_diox marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_dioxSource recognized as marginal(s).
## - - - Evaluating ERF_vit ...
##
## done(0 secs)!
## - - - Checking ERF_vit marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_vitSource recognized as marginal(s).
## - - - Evaluating ERF_micr ...
##
## done(0 secs)!
## - - - Checking ERF_micr marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_micrSource recognized as marginal(s).
## - - - Evaluating ERF_pfas ...
##
## done(0 secs)!
## - - - Checking ERF_pfas marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_pfasSource recognized as marginal(s).
## - - - Evaluating ERFchoice ...
##
## done(0 secs)!
## - - - Checking ERFchoice marginals ... Exposure_agent, Response, Scaling, Exposure, ER_function, ERFchoiceSource recognized as marginal(s).
##
## -- done(2.95 mins)!
## - - Checking ERF marginals ... Exposure_agent, Response, Exposure, ER_function, Scaling, Observation, ERFSource recognized as marginal(s).
## - - Processing ERF marginal collapses ... done!
## - - Evaluating RR ...
## - - - Processing dose marginal collapses ... done!
## - - - Processing ERF marginal collapses ... done!
##
## -- done(0.18 secs)!
## - - Checking RR marginals ... Exposure_agent, Response, ER_function, Scaling, ERFSource, Scenario, total_amountSource, Gender, fish_proportionSource, amountSource, conc_pfas_rawSource, doseSource, Age, RRSource recognized as marginal(s).
## - - Evaluating frexposed ...
##
## done(0 secs)!
## - - Checking frexposed marginals ... frexposedSource recognized as marginal(s).
## - - Evaluating P_illness ...
##
## done(0 secs)!
## - - Checking P_illness marginals ... Response, Illness, Age, P_illnessSource recognized as marginal(s).
##
## - done(8.1 mins)!
## - Checking PAF marginals ... Exposure_agent, Response, ER_function, Scaling, ERFSource, Scenario, total_amountSource, Gender, fish_proportionSource, amountSource, conc_pfas_rawSource, doseSource, frexposedSource, Age, incidenceSource, Adjust, RRSource, PAFSource recognized as marginal(s).
## - Processing PAF marginal collapses ...## Warning in oapply(variable, FUN = fun[[i]], cols = cols[[i]], na.rm = TRUE):
## While oapplying PAF, found NAs in indices: Age, Adjust. They were automatically
## filled using fillna, which may result in a multiplied population. Please check
## your ovariable before using oapply.## done!
##
## done(8.91 mins)!
## Checking BoDattr marginals ... Response, Age, Gender, Adjust, InpBoDSource, Exposure_agent, Scenario, total_amountSource, fish_proportionSource, amountSource, conc_pfas_rawSource, PAFSource, BoDattrSource recognized as marginal(s).oprint(summary(amount,"mean"))## Kala Scenario Gender Fish mean
## 1 Kaupallinen BAU Female Average fish 1.2697279
## 2 Muu tuonti BAU Female Average fish 11.1101191
## 3 Vapaa-ajan BAU Female Average fish 3.9282207
## 4 Kaupallinen BAU Male Average fish 1.9045919
## 5 Muu tuonti BAU Male Average fish 16.6651787
## 6 Vapaa-ajan BAU Male Average fish 5.8923310
## 7 Silakka BAU Female Herring 0.6348640
## 8 Silakka BAU Male Herring 0.9522959
## 9 Kirjolohi BAU Female Rainbow trout 2.6584928
## 10 Tuontikirjolohi BAU Female Rainbow trout 1.9045919
## 11 Kirjolohi BAU Male Rainbow trout 3.9877392
## 12 Tuontikirjolohi BAU Male Rainbow trout 2.8568878
## 13 Tuontilohi BAU Female Salmon 9.3642433
## 14 Tuontilohi BAU Male Salmon 14.0463649
## 15 Kasvatettu BAU Female Whitefish 0.2380740
## 16 Kasvatettu BAU Male Whitefish 0.3571110oprint(summary(BoD,marginals=c("Age","Response"),"mean"))## Age Response mean
## 1 0 - 4 Cancer morbidity 309.5200
## 2 10 - 14 Cancer morbidity 347.0750
## 3 15 - 19 Cancer morbidity 414.9750
## 4 20 - 24 Cancer morbidity 569.5250
## 5 25 - 29 Cancer morbidity 795.4550
## 6 30 - 34 Cancer morbidity 1147.9850
## 7 35 - 39 Cancer morbidity 1679.1250
## 8 40 - 44 Cancer morbidity 2452.5650
## 9 45 - 49 Cancer morbidity 3849.3800
## 10 5 - 9 Cancer morbidity 345.5200
## 11 50 - 54 Cancer morbidity 6987.9000
## 12 55 - 59 Cancer morbidity 11237.3300
## 13 60 - 64 Cancer morbidity 16331.9100
## 14 65 - 69 Cancer morbidity 21899.6500
## 15 70 - 74 Cancer morbidity 25551.8600
## 16 75 - 79 Cancer morbidity 18015.2300
## 17 80 - 84 Cancer morbidity 14650.9150
## 18 85 - 89 Cancer morbidity 8413.0550
## 19 90 - 94 Cancer morbidity 3537.1750
## 20 Undefined Dioxin recommendation tolerable daily intake 305.4067
## 21 Undefined Dioxin recommendation tolerable daily intake 2018 888.4560
## 22 0 - 4 Immunosuppression 295.9850
## 23 5 - 9 Immunosuppression 261.5350
## 24 0 - 4 Loss in child's IQ points 16778.3774
## 25 Undefined PFAS TWI 2776.4249
## 26 0 - 4 Sperm concentration 4478.6700
## 27 Undefined Vitamin D recommendation 610.8135
## 28 0 - 4 Yes or no dental defect 343.9619
## 29 0 - 4 All-cause mortality 4867.4650
## 30 10 - 14 All-cause mortality 1085.5450
## 31 15 - 19 All-cause mortality 2965.9450
## 32 20 - 24 All-cause mortality 5274.8450
## 33 25 - 29 All-cause mortality 6263.5400
## 34 30 - 34 All-cause mortality 6710.3100
## 35 35 - 39 All-cause mortality 8275.6450
## 36 40 - 44 All-cause mortality 10232.0150
## 37 45 - 49 All-cause mortality 13691.2800
## 38 5 - 9 All-cause mortality 931.1250
## 39 50 - 54 All-cause mortality 21321.8650
## 40 55 - 59 All-cause mortality 29861.9100
## 41 60 - 64 All-cause mortality 39200.8750
## 42 65 - 69 All-cause mortality 49295.5050
## 43 70 - 74 All-cause mortality 59799.7000
## 44 75 - 79 All-cause mortality 49428.5150
## 45 80 - 84 All-cause mortality 51968.9150
## 46 85 - 89 All-cause mortality 42765.3450
## 47 90 - 94 All-cause mortality 26756.6400
## 48 0 - 4 Depression 0.4150
## 49 10 - 14 Depression 555.8850
## 50 15 - 19 Depression 1320.8250
## 51 20 - 24 Depression 1841.0150
## 52 25 - 29 Depression 1801.4650
## 53 30 - 34 Depression 1590.0200
## 54 35 - 39 Depression 1721.7250
## 55 40 - 44 Depression 1700.2800
## 56 45 - 49 Depression 1534.1700
## 57 5 - 9 Depression 76.4850
## 58 50 - 54 Depression 1718.9550
## 59 55 - 59 Depression 1767.1650
## 60 60 - 64 Depression 1676.7400
## 61 65 - 69 Depression 1561.5750
## 62 70 - 74 Depression 1426.2350
## 63 75 - 79 Depression 799.7000
## 64 80 - 84 Depression 581.2950
## 65 85 - 89 Depression 363.0750
## 66 90 - 94 Depression 186.1150
## 67 0 - 4 CHD2 mortality 53.6100
## 68 10 - 14 CHD2 mortality 75.3150
## 69 15 - 19 CHD2 mortality 145.3450
## 70 20 - 24 CHD2 mortality 254.1000
## 71 25 - 29 CHD2 mortality 437.6450
## 72 30 - 34 CHD2 mortality 636.4200
## 73 35 - 39 CHD2 mortality 1168.2600
## 74 40 - 44 CHD2 mortality 1977.5650
## 75 45 - 49 CHD2 mortality 3117.3100
## 76 5 - 9 CHD2 mortality 48.7800
## 77 50 - 54 CHD2 mortality 5476.2900
## 78 55 - 59 CHD2 mortality 8747.4650
## 79 60 - 64 CHD2 mortality 13077.9700
## 80 65 - 69 CHD2 mortality 18100.7500
## 81 70 - 74 CHD2 mortality 24549.2850
## 82 75 - 79 CHD2 mortality 23413.4850
## 83 80 - 84 CHD2 mortality 27899.3350
## 84 85 - 89 CHD2 mortality 25240.0350
## 85 90 - 94 CHD2 mortality 16630.6150
## 86 15 - 19 Breast cancer 2.0900
## 87 20 - 24 Breast cancer 7.7850
## 88 25 - 29 Breast cancer 36.8200
## 89 30 - 34 Breast cancer 145.1600
## 90 35 - 39 Breast cancer 263.6000
## 91 40 - 44 Breast cancer 433.3500
## 92 45 - 49 Breast cancer 704.9550
## 93 50 - 54 Breast cancer 1045.3250
## 94 55 - 59 Breast cancer 1341.8150
## 95 60 - 64 Breast cancer 1513.2150
## 96 65 - 69 Breast cancer 1649.3300
## 97 70 - 74 Breast cancer 1678.0350
## 98 75 - 79 Breast cancer 1232.0100
## 99 80 - 84 Breast cancer 950.8800
## 100 85 - 89 Breast cancer 578.1900
## 101 90 - 94 Breast cancer 296.1800oprint(summary(BoDattr,marginals=c("Age","Response"),"mean"))## Age Response mean
## 1 0 - 4 Cancer morbidity 8.092797e+01
## 2 10 - 14 Cancer morbidity 8.850176e+01
## 3 15 - 19 Cancer morbidity 1.079079e+02
## 4 20 - 24 Cancer morbidity 1.483914e+02
## 5 25 - 29 Cancer morbidity 2.074377e+02
## 6 30 - 34 Cancer morbidity 2.931146e+02
## 7 35 - 39 Cancer morbidity 4.269816e+02
## 8 40 - 44 Cancer morbidity 6.183735e+02
## 9 45 - 49 Cancer morbidity 9.690704e+02
## 10 5 - 9 Cancer morbidity 8.825348e+01
## 11 50 - 54 Cancer morbidity 1.781467e+03
## 12 55 - 59 Cancer morbidity 2.926500e+03
## 13 60 - 64 Cancer morbidity 4.292083e+03
## 14 65 - 69 Cancer morbidity 5.788529e+03
## 15 70 - 74 Cancer morbidity 6.750697e+03
## 16 75 - 79 Cancer morbidity 4.726815e+03
## 17 80 - 84 Cancer morbidity 3.804571e+03
## 18 Undefined Dioxin recommendation tolerable daily intake 3.054067e+02
## 19 Undefined Dioxin recommendation tolerable daily intake 2018 8.884560e+02
## 20 0 - 4 Immunosuppression 6.818641e+00
## 21 5 - 9 Immunosuppression 7.294478e+00
## 22 0 - 4 Loss in child's IQ points 6.974660e+02
## 23 Undefined PFAS TWI 1.370935e+03
## 24 0 - 4 Sperm concentration 4.548773e+03
## 25 Undefined Vitamin D recommendation 6.108135e+02
## 26 0 - 4 Yes or no dental defect 2.530885e+03
## 27 0 - 4 All-cause mortality -3.940215e+02
## 28 10 - 14 All-cause mortality -8.729181e+01
## 29 15 - 19 All-cause mortality -2.510075e+02
## 30 20 - 24 All-cause mortality -4.582433e+02
## 31 25 - 29 All-cause mortality -5.472993e+02
## 32 30 - 34 All-cause mortality -5.805901e+02
## 33 35 - 39 All-cause mortality -7.164739e+02
## 34 40 - 44 All-cause mortality -8.793680e+02
## 35 45 - 49 All-cause mortality -1.166310e+03
## 36 5 - 9 All-cause mortality -7.424681e+01
## 37 50 - 54 All-cause mortality -1.808052e+03
## 38 55 - 59 All-cause mortality -2.535812e+03
## 39 60 - 64 All-cause mortality -3.319267e+03
## 40 65 - 69 All-cause mortality -4.145066e+03
## 41 70 - 74 All-cause mortality -4.980650e+03
## 42 75 - 79 All-cause mortality -4.044850e+03
## 43 80 - 84 All-cause mortality -4.134937e+03
## 44 0 - 4 Depression -7.647127e-02
## 45 10 - 14 Depression -9.937062e+01
## 46 15 - 19 Depression -2.347323e+02
## 47 20 - 24 Depression -3.276647e+02
## 48 25 - 29 Depression -3.226701e+02
## 49 30 - 34 Depression -2.865779e+02
## 50 35 - 39 Depression -3.110681e+02
## 51 40 - 44 Depression -3.073361e+02
## 52 45 - 49 Depression -2.766914e+02
## 53 5 - 9 Depression -1.387565e+01
## 54 50 - 54 Depression -3.090770e+02
## 55 55 - 59 Depression -3.166728e+02
## 56 60 - 64 Depression -2.997941e+02
## 57 65 - 69 Depression -2.791761e+02
## 58 70 - 74 Depression -2.551780e+02
## 59 75 - 79 Depression -1.425020e+02
## 60 80 - 84 Depression -1.024717e+02
## 61 0 - 4 CHD2 mortality -8.359787e+00
## 62 10 - 14 CHD2 mortality -1.170566e+01
## 63 15 - 19 CHD2 mortality -2.262663e+01
## 64 20 - 24 CHD2 mortality -3.966234e+01
## 65 25 - 29 CHD2 mortality -6.844631e+01
## 66 30 - 34 CHD2 mortality -9.964232e+01
## 67 35 - 39 CHD2 mortality -1.833859e+02
## 68 40 - 44 CHD2 mortality -3.106654e+02
## 69 45 - 49 CHD2 mortality -4.899542e+02
## 70 5 - 9 CHD2 mortality -7.579246e+00
## 71 50 - 54 CHD2 mortality -8.613391e+02
## 72 55 - 59 CHD2 mortality -1.375930e+03
## 73 60 - 64 CHD2 mortality -2.055932e+03
## 74 65 - 69 CHD2 mortality -2.841726e+03
## 75 70 - 74 CHD2 mortality -3.847268e+03
## 76 75 - 79 CHD2 mortality -3.660036e+03
## 77 80 - 84 CHD2 mortality -4.346514e+03
## 78 15 - 19 Breast cancer -4.624153e-01
## 79 20 - 24 Breast cancer -1.587543e+00
## 80 25 - 29 Breast cancer -7.326368e+00
## 81 30 - 34 Breast cancer -2.873340e+01
## 82 35 - 39 Breast cancer -5.218217e+01
## 83 40 - 44 Breast cancer -8.577757e+01
## 84 45 - 49 Breast cancer -1.395958e+02
## 85 50 - 54 Breast cancer -2.070689e+02
## 86 55 - 59 Breast cancer -2.658472e+02
## 87 60 - 64 Breast cancer -2.998331e+02
## 88 65 - 69 Breast cancer -3.267095e+02
## 89 70 - 74 Breast cancer -3.327787e+02
## 90 75 - 79 Breast cancer -2.441402e+02
## 91 80 - 84 Breast cancer -1.884558e+02oprint(summary(BoDattr,marginals=c("Exposure_agent","Response"),"mean"))## Exposure_agent Response
## 1 TEQ Cancer morbidity
## 2 TEQ Dioxin recommendation tolerable daily intake
## 3 TEQ Dioxin recommendation tolerable daily intake 2018
## 4 PFAS Immunosuppression
## 5 DHA Loss in child's IQ points
## 6 MeHg Loss in child's IQ points
## 7 PFAS PFAS TWI
## 8 TEQ Sperm concentration
## 9 Vitamin D Vitamin D recommendation
## 10 TEQ Yes or no dental defect
## 11 Fish All-cause mortality
## 12 Fish Depression
## 13 Omega3 CHD2 mortality
## 14 Omega3 Breast cancer
## mean
## 1 1947.036623
## 2 305.406736
## 3 888.455958
## 4 7.056559
## 5 -3292.013959
## 6 4686.945878
## 7 1370.934688
## 8 4548.773342
## 9 610.813471
## 10 2530.884510
## 11 -1771.969747
## 12 -228.525606
## 13 -1190.045495
## 14 -155.749901oprint(summary(BoDattr,marginals=c("Gender","Response"),"mean"))## Gender Response mean
## 1 Female Cancer morbidity 1378.857757
## 2 Male Cancer morbidity 2515.215489
## 3 Female Dioxin recommendation tolerable daily intake 309.207840
## 4 Male Dioxin recommendation tolerable daily intake 301.605631
## 5 Female Dioxin recommendation tolerable daily intake 2018 899.513715
## 6 Male Dioxin recommendation tolerable daily intake 2018 877.398200
## 7 Female Immunosuppression 5.205862
## 8 Male Immunosuppression 8.907256
## 9 Female Loss in child's IQ points 542.752772
## 10 Male Loss in child's IQ points 852.179148
## 11 Female PFAS TWI 0.000000
## 12 Male PFAS TWI 2741.869376
## 13 Female Sperm concentration 3539.756033
## 14 Male Sperm concentration 5557.790651
## 15 Female Vitamin D recommendation 618.415679
## 16 Male Vitamin D recommendation 603.211263
## 17 Female Yes or no dental defect 1969.479031
## 18 Male Yes or no dental defect 3092.289990
## 19 Female All-cause mortality -1003.248948
## 20 Male All-cause mortality -2540.690546
## 21 Female Depression -235.343301
## 22 Male Depression -221.707910
## 23 Female CHD2 mortality -815.643510
## 24 Male CHD2 mortality -1564.447480
## 25 Female Breast cancer -308.799252
## 26 Male Breast cancer -2.700550oprint(summary(case_burden,"mean"))## Response mean
## 1 Dioxin recommendation tolerable daily intake 0.001004988
## 2 Dioxin recommendation tolerable daily intake 2018 0.001004988
## 3 Loss in child's IQ points 0.110000000
## 4 PFAS TWI 0.001004988
## 5 Sperm concentration 2.500000000
## 6 Vitamin D recommendation 0.001004988
## 7 Yes or no dental defect 0.060000000oprint(summary(conc,"mean"))## Fish Compound mean
## 1 Average fish ALA 0.69000000
## 2 Bream ALA 0.22000000
## 3 Herring ALA 1.74000000
## 4 Pike ALA 0.08000000
## 5 Rainbow trout ALA 4.81000000
## 6 Roach ALA 0.10000000
## 7 Salmon ALA 7.96000000
## 8 Vendace ALA 1.35000000
## 9 Whitefish ALA 2.22000000
## 10 Average fish DHA 2.54000000
## 11 Bream DHA 2.73000000
## 12 Herring DHA 5.86000000
## 13 Pike DHA 0.30000000
## 14 Rainbow trout DHA 7.57000000
## 15 Roach DHA 2.87000000
## 16 Salmon DHA 6.69000000
## 17 Vendace DHA 3.00000000
## 18 Whitefish DHA 3.94000000
## 19 Average fish Fish 1.00000000
## 20 Bream Fish 1.00000000
## 21 Herring Fish 1.00000000
## 22 Pike Fish 1.00000000
## 23 Rainbow trout Fish 1.00000000
## 24 Roach Fish 1.00000000
## 25 Salmon Fish 1.00000000
## 26 Vendace Fish 1.00000000
## 27 Whitefish Fish 1.00000000
## 28 Average fish MeHg 0.11413395
## 29 Bream MeHg 0.63339772
## 30 Burbot MeHg 1.02686230
## 31 Chaetiliidae MeHg 0.01634649
## 32 Flounder MeHg 0.04582576
## 33 Fourhorn sculpin MeHg 0.64990137
## 34 Herring MeHg 0.02680994
## 35 Perch MeHg 0.14874312
## 36 Pike MeHg 0.90310953
## 37 Pike-perch MeHg 0.25432475
## 38 Roach MeHg 0.26739658
## 39 Sprat MeHg 0.03187976
## 40 White bream MeHg 1.21073890
## 41 Whitefish MeHg 0.06619502
## 42 Average fish Omega3 7.00000000
## 43 Bream Omega3 6.00000000
## 44 Herring Omega3 24.00000000
## 45 Pike Omega3 0.50000000
## 46 Rainbow trout Omega3 18.00000000
## 47 Roach Omega3 5.00000000
## 48 Salmon Omega3 23.00000000
## 49 Vendace Omega3 10.00000000
## 50 Whitefish Omega3 10.00000000
## 51 Average fish PFAS 8.79312500
## 52 Bream PFAS 4.95000000
## 53 Eel PFAS 7.89500000
## 54 Herring PFAS 1.65500000
## 55 Perch PFAS 8.79312500
## 56 Pike PFAS 8.79312500
## 57 Pike-perch PFAS 2.59500000
## 58 Rainbow trout PFAS 8.79312500
## 59 Roach PFAS 8.79312500
## 60 Salmon PFAS 8.79312500
## 61 Vendace PFAS 8.79312500
## 62 Whitefish PFAS 8.79312500
## 63 Baltic herring TEQ 90.08722156
## 64 Bream TEQ 6.20175433
## 65 Burbot TEQ 0.82549312
## 66 Cod TEQ 1.67367981
## 67 Flounder TEQ 15.36558024
## 68 Perch TEQ 3.56897420
## 69 Pike TEQ 2.14862962
## 70 Pike-perch TEQ 2.70814568
## 71 River lamprey TEQ 7.34301416
## 72 Roach TEQ 0.72124877
## 73 Salmon TEQ 20.00778296
## 74 Sea trout TEQ 9.31778459
## 75 Sprat TEQ 4.15456559
## 76 Vendace TEQ 0.99253888
## 77 Whitefish TEQ 5.74967872
## 78 Average fish Vitamin D 0.10500000
## 79 Bream Vitamin D 0.14000000
## 80 Herring Vitamin D 0.15600000
## 81 Pike Vitamin D 0.02100000
## 82 Rainbow trout Vitamin D 0.05100000
## 83 Roach Vitamin D 0.10000000
## 84 Salmon Vitamin D 0.06700000
## 85 Vendace Vitamin D 0.09400000
## 86 Whitefish Vitamin D 0.14400000oprint(summary(dose,"mean"))## Scenario Gender Exposure_agent Scaling mean
## 1 BAU Female ALA BW 1.56247954
## 2 BAU Male ALA BW 2.34371931
## 3 BAU Female DHA BW 2.04671635
## 4 BAU Male DHA BW 3.07007453
## 5 BAU Female Fish BW 0.44440477
## 6 BAU Male Fish BW 0.66660715
## 7 BAU Female MeHg BW 0.02705835
## 8 BAU Male MeHg BW 0.04058752
## 9 BAU Female Omega3 BW 6.13267239
## 10 BAU Male Omega3 BW 9.19900859
## 11 BAU Female PFAS BW 3.84296754
## 12 BAU Male PFAS BW 5.76445130
## 13 BAU Female TEQ BW 2.69609423
## 14 BAU Male TEQ BW 4.04414134
## 15 BAU Female Vitamin D BW 0.03865415
## 16 BAU Male Vitamin D BW 0.05798122
## 17 BAU Female ALA Log10 2.03891238
## 18 BAU Male ALA Log10 2.21500364
## 19 BAU Female DHA Log10 2.15615570
## 20 BAU Male DHA Log10 2.33224696
## 21 BAU Female Fish Log10 1.49287675
## 22 BAU Male Fish Log10 1.66896801
## 23 BAU Female MeHg Log10 0.27739929
## 24 BAU Male MeHg Log10 0.45349055
## 25 BAU Female Omega3 Log10 2.63274781
## 26 BAU Male Omega3 Log10 2.80883906
## 27 BAU Female PFAS Log10 2.42976476
## 28 BAU Male PFAS Log10 2.60585601
## 29 BAU Female TEQ Log10 2.27583311
## 30 BAU Male TEQ Log10 2.45192437
## 31 BAU Female Vitamin D Log10 0.43229411
## 32 BAU Male Vitamin D Log10 0.60838537
## 33 BAU Female ALA None 109.37356784
## 34 BAU Male ALA None 164.06035175
## 35 BAU Female DHA None 143.27014453
## 36 BAU Male DHA None 214.90521680
## 37 BAU Female Fish None 31.10833358
## 38 BAU Male Fish None 46.66250037
## 39 BAU Female MeHg None 1.89408422
## 40 BAU Male MeHg None 2.84112633
## 41 BAU Female Omega3 None 429.28706762
## 42 BAU Male Omega3 None 643.93060143
## 43 BAU Female PFAS None 269.00772748
## 44 BAU Male PFAS None 403.51159123
## 45 BAU Female TEQ None 188.72659592
## 46 BAU Male TEQ None 283.08989388
## 47 BAU Female Vitamin D None 2.70579016
## 48 BAU Male Vitamin D None 4.05868524oprint(summary(ERF,"mean"))## Exposure_agent Response
## 1 Fish All-cause mortality
## 2 Omega3 Breast cancer
## 3 Fish Depression
## 4 DHA Loss in child's IQ points
## 5 TEQ Sperm concentration
## 6 TEQ Yes or no dental defect
## 7 Omega3 CHD2 mortality
## 8 Vitamin D Vitamin D recommendation
## 9 MeHg Loss in child's IQ points
## 10 PFAS Immunosuppression
## 11 PFAS PFAS TWI
## 12 TEQ Cancer morbidity
## 13 TEQ Dioxin recommendation tolerable daily intake
## 14 TEQ Dioxin recommendation tolerable daily intake 2018
## 15 Fish All-cause mortality
## 16 Omega3 Breast cancer
## 17 Fish Depression
## 18 DHA Loss in child's IQ points
## 19 TEQ Sperm concentration
## 20 TEQ Yes or no dental defect
## 21 Omega3 CHD2 mortality
## 22 Vitamin D Vitamin D recommendation
## 23 MeHg Loss in child's IQ points
## 24 PFAS Immunosuppression
## 25 PFAS PFAS TWI
## 26 TEQ Cancer morbidity
## 27 TEQ Dioxin recommendation tolerable daily intake
## 28 TEQ Dioxin recommendation tolerable daily intake 2018
## ER_function Scaling Observation mean
## 1 RR None ERF 0.997871700
## 2 RR None ERF 0.999487200
## 3 RR None ERF 0.994690400
## 4 ERS None ERF -0.001300000
## 5 ERS None ERF 0.000060000
## 6 ERS None ERF 0.001390971
## 7 Relative Hill None ERF -0.170000000
## 8 Step None ERF 100.000000000
## 9 ERS BW ERF 9.800000000
## 10 ERS BW ERF 0.022700000
## 11 TWI BW ERF 4.400000000
## 12 CSF BW ERF 0.000500000
## 13 TDI BW ERF 2.000000000
## 14 TDI BW ERF 0.288900000
## 15 RR None Threshold 0.000000000
## 16 RR None Threshold 0.000000000
## 17 RR None Threshold 0.000000000
## 18 ERS None Threshold 0.000000000
## 19 ERS None Threshold 0.000000000
## 20 ERS None Threshold 0.000000000
## 21 Relative Hill None Threshold 47.000000000
## 22 Step None Threshold 10.000000000
## 23 ERS BW Threshold 0.000000000
## 24 ERS BW Threshold 0.000000000
## 25 TWI BW Threshold 0.000000000
## 26 CSF BW Threshold 0.000000000
## 27 TDI BW Threshold 0.000000000
## 28 TDI BW Threshold 0.000000000oprint(summary(expo_dir,"mean"))## Warning in oapply(conc, cols = "", FUN = mean): While oapplying conc, found NAs
## in indices: conc_pfas_rawSource. They were automatically filled using fillna,
## which may result in a multiplied population. Please check your ovariable before
## using oapply.## Fish Compound Kala Scenario Gender mean
## 1 Whitefish ALA Kasvatettu BAU Female 0.52852424
## 2 Whitefish DHA Kasvatettu BAU Female 0.93801149
## 3 Whitefish Fish Kasvatettu BAU Female 0.23807398
## 4 Whitefish MeHg Kasvatettu BAU Female 0.01575931
## 5 Whitefish Omega3 Kasvatettu BAU Female 2.38073981
## 6 Whitefish PFAS Kasvatettu BAU Female 2.09341428
## 7 Whitefish TEQ Kasvatettu BAU Female 1.36884891
## 8 Whitefish Vitamin D Kasvatettu BAU Female 0.03428265
## 9 Average fish ALA Kaupallinen BAU Female 0.87611225
## 10 Average fish DHA Kaupallinen BAU Female 3.22510887
## 11 Average fish Fish Kaupallinen BAU Female 1.26972790
## 12 Average fish MeHg Kaupallinen BAU Female 0.21890186
## 13 Average fish Omega3 Kaupallinen BAU Female 8.88809531
## 14 Average fish PFAS Kaupallinen BAU Female 8.91229950
## 15 Average fish TEQ Kaupallinen BAU Female 3.98511713
## 16 Average fish Vitamin D Kaupallinen BAU Female 0.13332143
## 17 Rainbow trout ALA Kirjolohi BAU Female 12.78735034
## 18 Rainbow trout DHA Kirjolohi BAU Female 20.12479045
## 19 Rainbow trout Fish Kirjolohi BAU Female 2.65849279
## 20 Rainbow trout Omega3 Kirjolohi BAU Female 47.85287028
## 21 Rainbow trout PFAS Kirjolohi BAU Female 23.37645944
## 22 Rainbow trout Vitamin D Kirjolohi BAU Female 0.13558313
## 23 Average fish ALA Muu tuonti BAU Female 7.66598220
## 24 Average fish DHA Muu tuonti BAU Female 28.21970261
## 25 Average fish Fish Muu tuonti BAU Female 11.11011914
## 26 Average fish MeHg Muu tuonti BAU Female 1.91539131
## 27 Average fish Omega3 Muu tuonti BAU Female 77.77083395
## 28 Average fish PFAS Muu tuonti BAU Female 77.98262060
## 29 Average fish TEQ Muu tuonti BAU Female 34.86977486
## 30 Average fish Vitamin D Muu tuonti BAU Female 1.16656251
## 31 Herring ALA Silakka BAU Female 1.10466327
## 32 Herring DHA Silakka BAU Female 3.72030275
## 33 Herring Fish Silakka BAU Female 0.63486395
## 34 Herring MeHg Silakka BAU Female 0.01702066
## 35 Herring Omega3 Silakka BAU Female 15.23673482
## 36 Herring PFAS Silakka BAU Female 1.05069984
## 37 Herring Vitamin D Silakka BAU Female 0.09903878
## 38 Rainbow trout ALA Tuontikirjolohi BAU Female 9.16108681
## 39 Rainbow trout DHA Tuontikirjolohi BAU Female 14.41776032
## 40 Rainbow trout Fish Tuontikirjolohi BAU Female 1.90459185
## 41 Rainbow trout Omega3 Tuontikirjolohi BAU Female 34.28265333
## 42 Rainbow trout PFAS Tuontikirjolohi BAU Female 16.74731423
## 43 Rainbow trout Vitamin D Tuontikirjolohi BAU Female 0.09713418
## 44 Salmon ALA Tuontilohi BAU Female 74.53937644
## 45 Salmon DHA Tuontilohi BAU Female 62.64678749
## 46 Salmon Fish Tuontilohi BAU Female 9.36424327
## 47 Salmon Omega3 Tuontilohi BAU Female 215.37759525
## 48 Salmon PFAS Tuontilohi BAU Female 82.34096162
## 49 Salmon TEQ Tuontilohi BAU Female 187.35774701
## 50 Salmon Vitamin D Tuontilohi BAU Female 0.62740430
## 51 Average fish ALA Vapaa-ajan BAU Female 2.71047228
## 52 Average fish DHA Vapaa-ajan BAU Female 9.97768056
## 53 Average fish Fish Vapaa-ajan BAU Female 3.92822069
## 54 Average fish MeHg Vapaa-ajan BAU Female 0.67722764
## 55 Average fish Omega3 Vapaa-ajan BAU Female 27.49754486
## 56 Average fish PFAS Vapaa-ajan BAU Female 27.57242657
## 57 Average fish TEQ Vapaa-ajan BAU Female 12.32895611
## 58 Average fish Vitamin D Vapaa-ajan BAU Female 0.41246317
## 59 Whitefish ALA Kasvatettu BAU Male 0.79278636
## 60 Whitefish DHA Kasvatettu BAU Male 1.40701723
## 61 Whitefish Fish Kasvatettu BAU Male 0.35711097
## 62 Whitefish MeHg Kasvatettu BAU Male 0.02363897
## 63 Whitefish Omega3 Kasvatettu BAU Male 3.57110972
## 64 Whitefish PFAS Kasvatettu BAU Male 3.14012142
## 65 Whitefish TEQ Kasvatettu BAU Male 2.05327336
## 66 Whitefish Vitamin D Kasvatettu BAU Male 0.05142398
## 67 Average fish ALA Kaupallinen BAU Male 1.31416838
## 68 Average fish DHA Kaupallinen BAU Male 4.83766330
## 69 Average fish Fish Kaupallinen BAU Male 1.90459185
## 70 Average fish MeHg Kaupallinen BAU Male 0.32835280
## 71 Average fish Omega3 Kaupallinen BAU Male 13.33214296
## 72 Average fish PFAS Kaupallinen BAU Male 13.36844925
## 73 Average fish TEQ Kaupallinen BAU Male 5.97767569
## 74 Average fish Vitamin D Kaupallinen BAU Male 0.19998214
## 75 Rainbow trout ALA Kirjolohi BAU Male 19.18102550
## 76 Rainbow trout DHA Kirjolohi BAU Male 30.18718567
## 77 Rainbow trout Fish Kirjolohi BAU Male 3.98773919
## 78 Rainbow trout Omega3 Kirjolohi BAU Male 71.77930542
## 79 Rainbow trout PFAS Kirjolohi BAU Male 35.06468916
## 80 Rainbow trout Vitamin D Kirjolohi BAU Male 0.20337470
## 81 Average fish ALA Muu tuonti BAU Male 11.49897331
## 82 Average fish DHA Muu tuonti BAU Male 42.32955391
## 83 Average fish Fish Muu tuonti BAU Male 16.66517870
## 84 Average fish MeHg Muu tuonti BAU Male 2.87308697
## 85 Average fish Omega3 Muu tuonti BAU Male 116.65625093
## 86 Average fish PFAS Muu tuonti BAU Male 116.97393090
## 87 Average fish TEQ Muu tuonti BAU Male 52.30466229
## 88 Average fish Vitamin D Muu tuonti BAU Male 1.74984376
## 89 Herring ALA Silakka BAU Male 1.65699491
## 90 Herring DHA Silakka BAU Male 5.58045413
## 91 Herring Fish Silakka BAU Male 0.95229593
## 92 Herring MeHg Silakka BAU Male 0.02553100
## 93 Herring Omega3 Silakka BAU Male 22.85510222
## 94 Herring PFAS Silakka BAU Male 1.57604976
## 95 Herring Vitamin D Silakka BAU Male 0.14855816
## 96 Rainbow trout ALA Tuontikirjolohi BAU Male 13.74163021
## 97 Rainbow trout DHA Tuontikirjolohi BAU Male 21.62664048
## 98 Rainbow trout Fish Tuontikirjolohi BAU Male 2.85688778
## 99 Rainbow trout Omega3 Tuontikirjolohi BAU Male 51.42398000
## 100 Rainbow trout PFAS Tuontikirjolohi BAU Male 25.12097134
## 101 Rainbow trout Vitamin D Tuontikirjolohi BAU Male 0.14570128
## 102 Salmon ALA Tuontilohi BAU Male 111.80906467
## 103 Salmon DHA Tuontilohi BAU Male 93.97018123
## 104 Salmon Fish Tuontilohi BAU Male 14.04636491
## 105 Salmon Omega3 Tuontilohi BAU Male 323.06639288
## 106 Salmon PFAS Tuontilohi BAU Male 123.51144243
## 107 Salmon TEQ Tuontilohi BAU Male 281.03662052
## 108 Salmon Vitamin D Tuontilohi BAU Male 0.94110645
## 109 Average fish ALA Vapaa-ajan BAU Male 4.06570842
## 110 Average fish DHA Vapaa-ajan BAU Male 14.96652085
## 111 Average fish Fish Vapaa-ajan BAU Male 5.89233104
## 112 Average fish MeHg Vapaa-ajan BAU Male 1.01584146
## 113 Average fish Omega3 Vapaa-ajan BAU Male 41.24631729
## 114 Average fish PFAS Vapaa-ajan BAU Male 41.35863985
## 115 Average fish TEQ Vapaa-ajan BAU Male 18.49343417
## 116 Average fish Vitamin D Vapaa-ajan BAU Male 0.61869476#oprint(summary(exposure,"mean"))
#oprint(summary(fish_proportion,"mean"))
oprint(summary(incidence,"mean"))## Response Age Adjust mean
## 1 Immunosuppression 0 - 4 BAU 4.8100
## 2 Loss in child's IQ points 0 - 4 BAU 1.1920
## 3 Sperm concentration 0 - 4 BAU 0.0140
## 4 Yes or no dental defect 0 - 4 BAU 0.0448
## 5 Immunosuppression 5 - 9 BAU 3.9500
## 6 Dioxin recommendation tolerable daily intake Undefined BAU 0.1100
## 7 Dioxin recommendation tolerable daily intake 2018 Undefined BAU 0.3200
## 8 PFAS TWI Undefined BAU 1.0000
## 9 Vitamin D recommendation Undefined BAU 0.2200oprint(summary(PAF,"mean"))## Exposure_agent Response Scenario
## 1 DHA Loss in child's IQ points BAU
## 2 MeHg Loss in child's IQ points BAU
## 3 TEQ Cancer morbidity BAU
## 4 TEQ Sperm concentration BAU
## 5 TEQ Yes or no dental defect BAU
## 6 PFAS Immunosuppression BAU
## 7 Fish All-cause mortality BAU
## 8 Omega3 Breast cancer BAU
## 9 Fish Depression BAU
## 10 Omega3 CHD2 mortality BAU
## 11 DHA Loss in child's IQ points BAU
## 12 MeHg Loss in child's IQ points BAU
## 13 TEQ Cancer morbidity BAU
## 14 TEQ Sperm concentration BAU
## 15 TEQ Yes or no dental defect BAU
## 16 PFAS Immunosuppression BAU
## 17 Fish All-cause mortality BAU
## 18 Omega3 Breast cancer BAU
## 19 Fish Depression BAU
## 20 Omega3 CHD2 mortality BAU
## 21 TEQ Cancer morbidity BAU
## 22 Fish All-cause mortality BAU
## 23 Omega3 Breast cancer BAU
## 24 Fish Depression BAU
## 25 Omega3 CHD2 mortality BAU
## 26 TEQ Cancer morbidity BAU
## 27 Fish All-cause mortality BAU
## 28 Omega3 Breast cancer BAU
## 29 Fish Depression BAU
## 30 Omega3 CHD2 mortality BAU
## 31 TEQ Cancer morbidity BAU
## 32 Fish All-cause mortality BAU
## 33 Omega3 Breast cancer BAU
## 34 Fish Depression BAU
## 35 Omega3 CHD2 mortality BAU
## 36 TEQ Cancer morbidity BAU
## 37 Fish All-cause mortality BAU
## 38 Omega3 Breast cancer BAU
## 39 Fish Depression BAU
## 40 Omega3 CHD2 mortality BAU
## 41 TEQ Cancer morbidity BAU
## 42 Fish All-cause mortality BAU
## 43 Omega3 Breast cancer BAU
## 44 Fish Depression BAU
## 45 Omega3 CHD2 mortality BAU
## 46 TEQ Cancer morbidity BAU
## 47 Fish All-cause mortality BAU
## 48 Omega3 Breast cancer BAU
## 49 Fish Depression BAU
## 50 Omega3 CHD2 mortality BAU
## 51 TEQ Cancer morbidity BAU
## 52 Fish All-cause mortality BAU
## 53 Omega3 Breast cancer BAU
## 54 Fish Depression BAU
## 55 Omega3 CHD2 mortality BAU
## 56 TEQ Cancer morbidity BAU
## 57 Fish All-cause mortality BAU
## 58 Omega3 Breast cancer BAU
## 59 Fish Depression BAU
## 60 Omega3 CHD2 mortality BAU
## 61 TEQ Cancer morbidity BAU
## 62 Fish All-cause mortality BAU
## 63 Omega3 Breast cancer BAU
## 64 Fish Depression BAU
## 65 Omega3 CHD2 mortality BAU
## 66 TEQ Cancer morbidity BAU
## 67 Fish All-cause mortality BAU
## 68 Omega3 Breast cancer BAU
## 69 Fish Depression BAU
## 70 Omega3 CHD2 mortality BAU
## 71 TEQ Cancer morbidity BAU
## 72 Fish All-cause mortality BAU
## 73 Omega3 Breast cancer BAU
## 74 Fish Depression BAU
## 75 Omega3 CHD2 mortality BAU
## 76 TEQ Cancer morbidity BAU
## 77 Fish All-cause mortality BAU
## 78 Omega3 Breast cancer BAU
## 79 Fish Depression BAU
## 80 Omega3 CHD2 mortality BAU
## 81 TEQ Cancer morbidity BAU
## 82 Fish All-cause mortality BAU
## 83 Omega3 Breast cancer BAU
## 84 Fish Depression BAU
## 85 Omega3 CHD2 mortality BAU
## 86 TEQ Cancer morbidity BAU
## 87 Fish All-cause mortality BAU
## 88 Omega3 Breast cancer BAU
## 89 Fish Depression BAU
## 90 Omega3 CHD2 mortality BAU
## 91 TEQ Cancer morbidity BAU
## 92 Fish All-cause mortality BAU
## 93 Omega3 Breast cancer BAU
## 94 Fish Depression BAU
## 95 Omega3 CHD2 mortality BAU
## 96 TEQ Cancer morbidity BAU
## 97 Fish All-cause mortality BAU
## 98 Omega3 Breast cancer BAU
## 99 Fish Depression BAU
## 100 Omega3 CHD2 mortality BAU
## 101 TEQ Cancer morbidity BAU
## 102 PFAS Immunosuppression BAU
## 103 Fish All-cause mortality BAU
## 104 Omega3 Breast cancer BAU
## 105 Fish Depression BAU
## 106 Omega3 CHD2 mortality BAU
## 107 TEQ Cancer morbidity BAU
## 108 PFAS Immunosuppression BAU
## 109 Fish All-cause mortality BAU
## 110 Omega3 Breast cancer BAU
## 111 Fish Depression BAU
## 112 Omega3 CHD2 mortality BAU
## 113 TEQ Cancer morbidity BAU
## 114 Fish All-cause mortality BAU
## 115 Omega3 Breast cancer BAU
## 116 Fish Depression BAU
## 117 Omega3 CHD2 mortality BAU
## 118 TEQ Cancer morbidity BAU
## 119 Fish All-cause mortality BAU
## 120 Omega3 Breast cancer BAU
## 121 Fish Depression BAU
## 122 Omega3 CHD2 mortality BAU
## 123 TEQ Cancer morbidity BAU
## 124 Fish All-cause mortality BAU
## 125 Omega3 Breast cancer BAU
## 126 Fish Depression BAU
## 127 Omega3 CHD2 mortality BAU
## 128 TEQ Cancer morbidity BAU
## 129 Fish All-cause mortality BAU
## 130 Omega3 Breast cancer BAU
## 131 Fish Depression BAU
## 132 Omega3 CHD2 mortality BAU
## 133 TEQ Cancer morbidity BAU
## 134 Fish All-cause mortality BAU
## 135 Omega3 Breast cancer BAU
## 136 Fish Depression BAU
## 137 Omega3 CHD2 mortality BAU
## 138 TEQ Cancer morbidity BAU
## 139 Fish All-cause mortality BAU
## 140 Omega3 Breast cancer BAU
## 141 Fish Depression BAU
## 142 Omega3 CHD2 mortality BAU
## 143 TEQ Cancer morbidity BAU
## 144 Fish All-cause mortality BAU
## 145 Omega3 Breast cancer BAU
## 146 Fish Depression BAU
## 147 Omega3 CHD2 mortality BAU
## 148 TEQ Cancer morbidity BAU
## 149 Fish All-cause mortality BAU
## 150 Omega3 Breast cancer BAU
## 151 Fish Depression BAU
## 152 Omega3 CHD2 mortality BAU
## 153 TEQ Cancer morbidity BAU
## 154 Fish All-cause mortality BAU
## 155 Omega3 Breast cancer BAU
## 156 Fish Depression BAU
## 157 Omega3 CHD2 mortality BAU
## 158 TEQ Cancer morbidity BAU
## 159 Fish All-cause mortality BAU
## 160 Omega3 Breast cancer BAU
## 161 Fish Depression BAU
## 162 Omega3 CHD2 mortality BAU
## 163 TEQ Cancer morbidity BAU
## 164 Fish All-cause mortality BAU
## 165 Omega3 Breast cancer BAU
## 166 Fish Depression BAU
## 167 Omega3 CHD2 mortality BAU
## 168 TEQ Cancer morbidity BAU
## 169 Fish All-cause mortality BAU
## 170 Omega3 Breast cancer BAU
## 171 Fish Depression BAU
## 172 Omega3 CHD2 mortality BAU
## 173 TEQ Cancer morbidity BAU
## 174 Fish All-cause mortality BAU
## 175 Omega3 Breast cancer BAU
## 176 Fish Depression BAU
## 177 Omega3 CHD2 mortality BAU
## 178 TEQ Cancer morbidity BAU
## 179 Fish All-cause mortality BAU
## 180 Omega3 Breast cancer BAU
## 181 Fish Depression BAU
## 182 Omega3 CHD2 mortality BAU
## 183 TEQ Cancer morbidity BAU
## 184 Fish All-cause mortality BAU
## 185 Omega3 Breast cancer BAU
## 186 Fish Depression BAU
## 187 Omega3 CHD2 mortality BAU
## 188 TEQ Cancer morbidity BAU
## 189 Fish All-cause mortality BAU
## 190 Omega3 Breast cancer BAU
## 191 Fish Depression BAU
## 192 Omega3 CHD2 mortality BAU
## 193 TEQ Cancer morbidity BAU
## 194 TEQ Dioxin recommendation tolerable daily intake BAU
## 195 TEQ Dioxin recommendation tolerable daily intake 2018 BAU
## 196 Vitamin D Vitamin D recommendation BAU
## 197 PFAS PFAS TWI BAU
## 198 Fish All-cause mortality BAU
## 199 Omega3 Breast cancer BAU
## 200 Fish Depression BAU
## 201 Omega3 CHD2 mortality BAU
## 202 TEQ Cancer morbidity BAU
## 203 TEQ Dioxin recommendation tolerable daily intake BAU
## 204 TEQ Dioxin recommendation tolerable daily intake 2018 BAU
## 205 Vitamin D Vitamin D recommendation BAU
## 206 PFAS PFAS TWI BAU
## 207 Fish All-cause mortality BAU
## 208 Omega3 Breast cancer BAU
## 209 Fish Depression BAU
## 210 Omega3 CHD2 mortality BAU
## Gender Age Adjust mean
## 1 Female 0 - 4 BAU -0.15625100
## 2 Female 0 - 4 BAU 0.22245956
## 3 Female 0 - 4 BAU 0.20518221
## 4 Female 0 - 4 BAU 0.80882827
## 5 Female 0 - 4 BAU 5.85966964
## 6 Female 0 - 4 BAU 0.01813625
## 7 Female 0 - 4 BAU -0.06412974
## 8 Female 0 - 4 BAU -0.19763757
## 9 Female 0 - 4 BAU -0.15262578
## 10 Female 0 - 4 BAU -0.15322440
## 11 Male 0 - 4 BAU -0.23437649
## 12 Male 0 - 4 BAU 0.33368933
## 13 Male 0 - 4 BAU 0.30777331
## 14 Male 0 - 4 BAU 1.21324240
## 15 Male 0 - 4 BAU 8.78950446
## 16 Male 0 - 4 BAU 0.02720438
## 17 Male 0 - 4 BAU -0.09463548
## 18 Male 0 - 4 BAU -0.28128638
## 19 Male 0 - 4 BAU -0.21996719
## 20 Male 0 - 4 BAU -0.15843589
## 21 Female 10 - 14 BAU 0.20518221
## 22 Female 10 - 14 BAU -0.06412974
## 23 Female 10 - 14 BAU -0.19763757
## 24 Female 10 - 14 BAU -0.15262578
## 25 Female 10 - 14 BAU -0.15322440
## 26 Male 10 - 14 BAU 0.30777331
## 27 Male 10 - 14 BAU -0.09463548
## 28 Male 10 - 14 BAU -0.28128638
## 29 Male 10 - 14 BAU -0.21996719
## 30 Male 10 - 14 BAU -0.15843589
## 31 Female 15 - 19 BAU 0.20518221
## 32 Female 15 - 19 BAU -0.06412974
## 33 Female 15 - 19 BAU -0.19763757
## 34 Female 15 - 19 BAU -0.15262578
## 35 Female 15 - 19 BAU -0.15322440
## 36 Male 15 - 19 BAU 0.30777331
## 37 Male 15 - 19 BAU -0.09463548
## 38 Male 15 - 19 BAU -0.28128638
## 39 Male 15 - 19 BAU -0.21996719
## 40 Male 15 - 19 BAU -0.15843589
## 41 Female 20 - 24 BAU 0.20518221
## 42 Female 20 - 24 BAU -0.06412974
## 43 Female 20 - 24 BAU -0.19763757
## 44 Female 20 - 24 BAU -0.15262578
## 45 Female 20 - 24 BAU -0.15322440
## 46 Male 20 - 24 BAU 0.30777331
## 47 Male 20 - 24 BAU -0.09463548
## 48 Male 20 - 24 BAU -0.28128638
## 49 Male 20 - 24 BAU -0.21996719
## 50 Male 20 - 24 BAU -0.15843589
## 51 Female 25 - 29 BAU 0.20518221
## 52 Female 25 - 29 BAU -0.06412974
## 53 Female 25 - 29 BAU -0.19763757
## 54 Female 25 - 29 BAU -0.15262578
## 55 Female 25 - 29 BAU -0.15322440
## 56 Male 25 - 29 BAU 0.30777331
## 57 Male 25 - 29 BAU -0.09463548
## 58 Male 25 - 29 BAU -0.28128638
## 59 Male 25 - 29 BAU -0.21996719
## 60 Male 25 - 29 BAU -0.15843589
## 61 Female 30 - 34 BAU 0.20518221
## 62 Female 30 - 34 BAU -0.06412974
## 63 Female 30 - 34 BAU -0.19763757
## 64 Female 30 - 34 BAU -0.15262578
## 65 Female 30 - 34 BAU -0.15322440
## 66 Male 30 - 34 BAU 0.30777331
## 67 Male 30 - 34 BAU -0.09463548
## 68 Male 30 - 34 BAU -0.28128638
## 69 Male 30 - 34 BAU -0.21996719
## 70 Male 30 - 34 BAU -0.15843589
## 71 Female 35 - 39 BAU 0.20518221
## 72 Female 35 - 39 BAU -0.06412974
## 73 Female 35 - 39 BAU -0.19763757
## 74 Female 35 - 39 BAU -0.15262578
## 75 Female 35 - 39 BAU -0.15322440
## 76 Male 35 - 39 BAU 0.30777331
## 77 Male 35 - 39 BAU -0.09463548
## 78 Male 35 - 39 BAU -0.28128638
## 79 Male 35 - 39 BAU -0.21996719
## 80 Male 35 - 39 BAU -0.15843589
## 81 Female 40 - 44 BAU 0.20518221
## 82 Female 40 - 44 BAU -0.06412974
## 83 Female 40 - 44 BAU -0.19763757
## 84 Female 40 - 44 BAU -0.15262578
## 85 Female 40 - 44 BAU -0.15322440
## 86 Male 40 - 44 BAU 0.30777331
## 87 Male 40 - 44 BAU -0.09463548
## 88 Male 40 - 44 BAU -0.28128638
## 89 Male 40 - 44 BAU -0.21996719
## 90 Male 40 - 44 BAU -0.15843589
## 91 Female 45 - 49 BAU 0.20518221
## 92 Female 45 - 49 BAU -0.06412974
## 93 Female 45 - 49 BAU -0.19763757
## 94 Female 45 - 49 BAU -0.15262578
## 95 Female 45 - 49 BAU -0.15322440
## 96 Male 45 - 49 BAU 0.30777331
## 97 Male 45 - 49 BAU -0.09463548
## 98 Male 45 - 49 BAU -0.28128638
## 99 Male 45 - 49 BAU -0.21996719
## 100 Male 45 - 49 BAU -0.15843589
## 101 Female 5 - 9 BAU 0.20518221
## 102 Female 5 - 9 BAU 0.02208490
## 103 Female 5 - 9 BAU -0.06412974
## 104 Female 5 - 9 BAU -0.19763757
## 105 Female 5 - 9 BAU -0.15262578
## 106 Female 5 - 9 BAU -0.15322440
## 107 Male 5 - 9 BAU 0.30777331
## 108 Male 5 - 9 BAU 0.03312735
## 109 Male 5 - 9 BAU -0.09463548
## 110 Male 5 - 9 BAU -0.28128638
## 111 Male 5 - 9 BAU -0.21996719
## 112 Male 5 - 9 BAU -0.15843589
## 113 Female 50 - 54 BAU 0.20518221
## 114 Female 50 - 54 BAU -0.06412974
## 115 Female 50 - 54 BAU -0.19763757
## 116 Female 50 - 54 BAU -0.15262578
## 117 Female 50 - 54 BAU -0.15322440
## 118 Male 50 - 54 BAU 0.30777331
## 119 Male 50 - 54 BAU -0.09463548
## 120 Male 50 - 54 BAU -0.28128638
## 121 Male 50 - 54 BAU -0.21996719
## 122 Male 50 - 54 BAU -0.15843589
## 123 Female 55 - 59 BAU 0.20518221
## 124 Female 55 - 59 BAU -0.06412974
## 125 Female 55 - 59 BAU -0.19763757
## 126 Female 55 - 59 BAU -0.15262578
## 127 Female 55 - 59 BAU -0.15322440
## 128 Male 55 - 59 BAU 0.30777331
## 129 Male 55 - 59 BAU -0.09463548
## 130 Male 55 - 59 BAU -0.28128638
## 131 Male 55 - 59 BAU -0.21996719
## 132 Male 55 - 59 BAU -0.15843589
## 133 Female 60 - 64 BAU 0.20518221
## 134 Female 60 - 64 BAU -0.06412974
## 135 Female 60 - 64 BAU -0.19763757
## 136 Female 60 - 64 BAU -0.15262578
## 137 Female 60 - 64 BAU -0.15322440
## 138 Male 60 - 64 BAU 0.30777331
## 139 Male 60 - 64 BAU -0.09463548
## 140 Male 60 - 64 BAU -0.28128638
## 141 Male 60 - 64 BAU -0.21996719
## 142 Male 60 - 64 BAU -0.15843589
## 143 Female 65 - 69 BAU 0.20518221
## 144 Female 65 - 69 BAU -0.06412974
## 145 Female 65 - 69 BAU -0.19763757
## 146 Female 65 - 69 BAU -0.15262578
## 147 Female 65 - 69 BAU -0.15322440
## 148 Male 65 - 69 BAU 0.30777331
## 149 Male 65 - 69 BAU -0.09463548
## 150 Male 65 - 69 BAU -0.28128638
## 151 Male 65 - 69 BAU -0.21996719
## 152 Male 65 - 69 BAU -0.15843589
## 153 Female 70 - 74 BAU 0.20518221
## 154 Female 70 - 74 BAU -0.06412974
## 155 Female 70 - 74 BAU -0.19763757
## 156 Female 70 - 74 BAU -0.15262578
## 157 Female 70 - 74 BAU -0.15322440
## 158 Male 70 - 74 BAU 0.30777331
## 159 Male 70 - 74 BAU -0.09463548
## 160 Male 70 - 74 BAU -0.28128638
## 161 Male 70 - 74 BAU -0.21996719
## 162 Male 70 - 74 BAU -0.15843589
## 163 Female 75 - 79 BAU 0.20518221
## 164 Female 75 - 79 BAU -0.06412974
## 165 Female 75 - 79 BAU -0.19763757
## 166 Female 75 - 79 BAU -0.15262578
## 167 Female 75 - 79 BAU -0.15322440
## 168 Male 75 - 79 BAU 0.30777331
## 169 Male 75 - 79 BAU -0.09463548
## 170 Male 75 - 79 BAU -0.28128638
## 171 Male 75 - 79 BAU -0.21996719
## 172 Male 75 - 79 BAU -0.15843589
## 173 Female 80 - 84 BAU 0.20518221
## 174 Female 80 - 84 BAU -0.06412974
## 175 Female 80 - 84 BAU -0.19763757
## 176 Female 80 - 84 BAU -0.15262578
## 177 Female 80 - 84 BAU -0.15322440
## 178 Male 80 - 84 BAU 0.30777331
## 179 Male 80 - 84 BAU -0.09463548
## 180 Male 80 - 84 BAU -0.28128638
## 181 Male 80 - 84 BAU -0.21996719
## 182 Male 80 - 84 BAU -0.15843589
## 183 Female 85+ BAU 0.20518221
## 184 Female 85+ BAU -0.06412974
## 185 Female 85+ BAU -0.19763757
## 186 Female 85+ BAU -0.15262578
## 187 Female 85+ BAU -0.15322440
## 188 Male 85+ BAU 0.30777331
## 189 Male 85+ BAU -0.09463548
## 190 Male 85+ BAU -0.28128638
## 191 Male 85+ BAU -0.21996719
## 192 Male 85+ BAU -0.15843589
## 193 Female Undefined BAU 0.20518221
## 194 Female Undefined BAU 1.00000000
## 195 Female Undefined BAU 1.00000000
## 196 Female Undefined BAU 1.00000000
## 197 Female Undefined BAU 0.00000000
## 198 Female Undefined BAU -0.06412974
## 199 Female Undefined BAU -0.19763757
## 200 Female Undefined BAU -0.15262578
## 201 Female Undefined BAU -0.15322440
## 202 Male Undefined BAU 0.30777331
## 203 Male Undefined BAU 1.00000000
## 204 Male Undefined BAU 1.00000000
## 205 Male Undefined BAU 1.00000000
## 206 Male Undefined BAU 1.00000000
## 207 Male Undefined BAU -0.09463548
## 208 Male Undefined BAU -0.28128638
## 209 Male Undefined BAU -0.21996719
## 210 Male Undefined BAU -0.15843589oprint(summary(population,"mean"))## Gender Age mean
## 1 Female 0 - 4 125040
## 2 Male 0 - 4 130884
## 3 Female 5 - 9 149633
## 4 Male 5 - 9 156654
## 5 Female 10 - 14 151113
## 6 Male 10 - 14 157712
## 7 Female 15 - 19 144441
## 8 Male 15 - 19 152230
## 9 Female 20 - 24 152265
## 10 Male 20 - 24 161679
## 11 Female 25 - 29 172593
## 12 Male 25 - 29 183092
## 13 Female 30 - 34 169653
## 14 Male 30 - 34 181115
## 15 Female 35 - 39 174660
## 16 Male 35 - 39 186122
## 17 Female 40 - 44 168547
## 18 Male 40 - 44 177928
## 19 Female 45 - 49 154391
## 20 Male 45 - 49 159982
## 21 Female 50 - 54 176612
## 22 Male 50 - 54 179182
## 23 Female 55 - 59 185152
## 24 Male 55 - 59 183719
## 25 Female 60 - 64 183336
## 26 Male 60 - 64 176283
## 27 Female 65 - 69 185685
## 28 Male 65 - 69 171275
## 29 Female 70 - 74 186034
## 30 Male 70 - 74 163697
## 31 Female 75 - 79 118190
## 32 Male 75 - 79 93987
## 33 Female 80 - 84 96256
## 34 Male 80 - 84 65140
## 35 Female 85+ 103429
## 36 Male 85+ 47581
## 37 Female Undefined 2797030
## 38 Male Undefined 2728262oprint(summary(RR,"mean"))## Exposure_agent Response ER_function Scaling Scenario Gender
## 1 Fish All-cause mortality RR None BAU Female
## 2 Omega3 Breast cancer RR None BAU Female
## 3 Fish Depression RR None BAU Female
## 4 Omega3 CHD2 mortality Relative Hill None BAU Female
## 5 Fish All-cause mortality RR None BAU Male
## 6 Omega3 Breast cancer RR None BAU Male
## 7 Fish Depression RR None BAU Male
## 8 Omega3 CHD2 mortality Relative Hill None BAU Male
## 9 Fish All-cause mortality RR None BAU Female
## 10 Omega3 Breast cancer RR None BAU Female
## 11 Fish Depression RR None BAU Female
## 12 Omega3 CHD2 mortality Relative Hill None BAU Female
## 13 Fish All-cause mortality RR None BAU Male
## 14 Omega3 Breast cancer RR None BAU Male
## 15 Fish Depression RR None BAU Male
## 16 Omega3 CHD2 mortality Relative Hill None BAU Male
## 17 Fish All-cause mortality RR None BAU Female
## 18 Omega3 Breast cancer RR None BAU Female
## 19 Fish Depression RR None BAU Female
## 20 Omega3 CHD2 mortality Relative Hill None BAU Female
## 21 Fish All-cause mortality RR None BAU Male
## 22 Omega3 Breast cancer RR None BAU Male
## 23 Fish Depression RR None BAU Male
## 24 Omega3 CHD2 mortality Relative Hill None BAU Male
## 25 Fish All-cause mortality RR None BAU Female
## 26 Omega3 Breast cancer RR None BAU Female
## 27 Fish Depression RR None BAU Female
## 28 Omega3 CHD2 mortality Relative Hill None BAU Female
## 29 Fish All-cause mortality RR None BAU Male
## 30 Omega3 Breast cancer RR None BAU Male
## 31 Fish Depression RR None BAU Male
## 32 Omega3 CHD2 mortality Relative Hill None BAU Male
## 33 Fish All-cause mortality RR None BAU Female
## 34 Omega3 Breast cancer RR None BAU Female
## 35 Fish Depression RR None BAU Female
## 36 Omega3 CHD2 mortality Relative Hill None BAU Female
## 37 Fish All-cause mortality RR None BAU Male
## 38 Omega3 Breast cancer RR None BAU Male
## 39 Fish Depression RR None BAU Male
## 40 Omega3 CHD2 mortality Relative Hill None BAU Male
## 41 Fish All-cause mortality RR None BAU Female
## 42 Omega3 Breast cancer RR None BAU Female
## 43 Fish Depression RR None BAU Female
## 44 Omega3 CHD2 mortality Relative Hill None BAU Female
## 45 Fish All-cause mortality RR None BAU Male
## 46 Omega3 Breast cancer RR None BAU Male
## 47 Fish Depression RR None BAU Male
## 48 Omega3 CHD2 mortality Relative Hill None BAU Male
## 49 Fish All-cause mortality RR None BAU Female
## 50 Omega3 Breast cancer RR None BAU Female
## 51 Fish Depression RR None BAU Female
## 52 Omega3 CHD2 mortality Relative Hill None BAU Female
## 53 Fish All-cause mortality RR None BAU Male
## 54 Omega3 Breast cancer RR None BAU Male
## 55 Fish Depression RR None BAU Male
## 56 Omega3 CHD2 mortality Relative Hill None BAU Male
## 57 Fish All-cause mortality RR None BAU Female
## 58 Omega3 Breast cancer RR None BAU Female
## 59 Fish Depression RR None BAU Female
## 60 Omega3 CHD2 mortality Relative Hill None BAU Female
## 61 Fish All-cause mortality RR None BAU Male
## 62 Omega3 Breast cancer RR None BAU Male
## 63 Fish Depression RR None BAU Male
## 64 Omega3 CHD2 mortality Relative Hill None BAU Male
## 65 Fish All-cause mortality RR None BAU Female
## 66 Omega3 Breast cancer RR None BAU Female
## 67 Fish Depression RR None BAU Female
## 68 Omega3 CHD2 mortality Relative Hill None BAU Female
## 69 Fish All-cause mortality RR None BAU Male
## 70 Omega3 Breast cancer RR None BAU Male
## 71 Fish Depression RR None BAU Male
## 72 Omega3 CHD2 mortality Relative Hill None BAU Male
## 73 Fish All-cause mortality RR None BAU Female
## 74 Omega3 Breast cancer RR None BAU Female
## 75 Fish Depression RR None BAU Female
## 76 Omega3 CHD2 mortality Relative Hill None BAU Female
## 77 Fish All-cause mortality RR None BAU Male
## 78 Omega3 Breast cancer RR None BAU Male
## 79 Fish Depression RR None BAU Male
## 80 Omega3 CHD2 mortality Relative Hill None BAU Male
## 81 Fish All-cause mortality RR None BAU Female
## 82 Omega3 Breast cancer RR None BAU Female
## 83 Fish Depression RR None BAU Female
## 84 Omega3 CHD2 mortality Relative Hill None BAU Female
## 85 Fish All-cause mortality RR None BAU Male
## 86 Omega3 Breast cancer RR None BAU Male
## 87 Fish Depression RR None BAU Male
## 88 Omega3 CHD2 mortality Relative Hill None BAU Male
## 89 Fish All-cause mortality RR None BAU Female
## 90 Omega3 Breast cancer RR None BAU Female
## 91 Fish Depression RR None BAU Female
## 92 Omega3 CHD2 mortality Relative Hill None BAU Female
## 93 Fish All-cause mortality RR None BAU Male
## 94 Omega3 Breast cancer RR None BAU Male
## 95 Fish Depression RR None BAU Male
## 96 Omega3 CHD2 mortality Relative Hill None BAU Male
## 97 Fish All-cause mortality RR None BAU Female
## 98 Omega3 Breast cancer RR None BAU Female
## 99 Fish Depression RR None BAU Female
## 100 Omega3 CHD2 mortality Relative Hill None BAU Female
## 101 Fish All-cause mortality RR None BAU Male
## 102 Omega3 Breast cancer RR None BAU Male
## 103 Fish Depression RR None BAU Male
## 104 Omega3 CHD2 mortality Relative Hill None BAU Male
## 105 Fish All-cause mortality RR None BAU Female
## 106 Omega3 Breast cancer RR None BAU Female
## 107 Fish Depression RR None BAU Female
## 108 Omega3 CHD2 mortality Relative Hill None BAU Female
## 109 Fish All-cause mortality RR None BAU Male
## 110 Omega3 Breast cancer RR None BAU Male
## 111 Fish Depression RR None BAU Male
## 112 Omega3 CHD2 mortality Relative Hill None BAU Male
## 113 Fish All-cause mortality RR None BAU Female
## 114 Omega3 Breast cancer RR None BAU Female
## 115 Fish Depression RR None BAU Female
## 116 Omega3 CHD2 mortality Relative Hill None BAU Female
## 117 Fish All-cause mortality RR None BAU Male
## 118 Omega3 Breast cancer RR None BAU Male
## 119 Fish Depression RR None BAU Male
## 120 Omega3 CHD2 mortality Relative Hill None BAU Male
## 121 Fish All-cause mortality RR None BAU Female
## 122 Omega3 Breast cancer RR None BAU Female
## 123 Fish Depression RR None BAU Female
## 124 Omega3 CHD2 mortality Relative Hill None BAU Female
## 125 Fish All-cause mortality RR None BAU Male
## 126 Omega3 Breast cancer RR None BAU Male
## 127 Fish Depression RR None BAU Male
## 128 Omega3 CHD2 mortality Relative Hill None BAU Male
## 129 Fish All-cause mortality RR None BAU Female
## 130 Omega3 Breast cancer RR None BAU Female
## 131 Fish Depression RR None BAU Female
## 132 Omega3 CHD2 mortality Relative Hill None BAU Female
## 133 Fish All-cause mortality RR None BAU Male
## 134 Omega3 Breast cancer RR None BAU Male
## 135 Fish Depression RR None BAU Male
## 136 Omega3 CHD2 mortality Relative Hill None BAU Male
## 137 Fish All-cause mortality RR None BAU Female
## 138 Omega3 Breast cancer RR None BAU Female
## 139 Fish Depression RR None BAU Female
## 140 Omega3 CHD2 mortality Relative Hill None BAU Female
## 141 Fish All-cause mortality RR None BAU Male
## 142 Omega3 Breast cancer RR None BAU Male
## 143 Fish Depression RR None BAU Male
## 144 Omega3 CHD2 mortality Relative Hill None BAU Male
## 145 Fish All-cause mortality RR None BAU Female
## 146 Omega3 Breast cancer RR None BAU Female
## 147 Fish Depression RR None BAU Female
## 148 Omega3 CHD2 mortality Relative Hill None BAU Female
## 149 Fish All-cause mortality RR None BAU Male
## 150 Omega3 Breast cancer RR None BAU Male
## 151 Fish Depression RR None BAU Male
## 152 Omega3 CHD2 mortality Relative Hill None BAU Male
## Age mean
## 1 0 - 4 0.9358703
## 2 0 - 4 0.8023624
## 3 0 - 4 0.8473742
## 4 0 - 4 0.8467756
## 5 0 - 4 0.9053645
## 6 0 - 4 0.7187136
## 7 0 - 4 0.7800328
## 8 0 - 4 0.8415641
## 9 10 - 14 0.9358703
## 10 10 - 14 0.8023624
## 11 10 - 14 0.8473742
## 12 10 - 14 0.8467756
## 13 10 - 14 0.9053645
## 14 10 - 14 0.7187136
## 15 10 - 14 0.7800328
## 16 10 - 14 0.8415641
## 17 15 - 19 0.9358703
## 18 15 - 19 0.8023624
## 19 15 - 19 0.8473742
## 20 15 - 19 0.8467756
## 21 15 - 19 0.9053645
## 22 15 - 19 0.7187136
## 23 15 - 19 0.7800328
## 24 15 - 19 0.8415641
## 25 20 - 24 0.9358703
## 26 20 - 24 0.8023624
## 27 20 - 24 0.8473742
## 28 20 - 24 0.8467756
## 29 20 - 24 0.9053645
## 30 20 - 24 0.7187136
## 31 20 - 24 0.7800328
## 32 20 - 24 0.8415641
## 33 25 - 29 0.9358703
## 34 25 - 29 0.8023624
## 35 25 - 29 0.8473742
## 36 25 - 29 0.8467756
## 37 25 - 29 0.9053645
## 38 25 - 29 0.7187136
## 39 25 - 29 0.7800328
## 40 25 - 29 0.8415641
## 41 30 - 34 0.9358703
## 42 30 - 34 0.8023624
## 43 30 - 34 0.8473742
## 44 30 - 34 0.8467756
## 45 30 - 34 0.9053645
## 46 30 - 34 0.7187136
## 47 30 - 34 0.7800328
## 48 30 - 34 0.8415641
## 49 35 - 39 0.9358703
## 50 35 - 39 0.8023624
## 51 35 - 39 0.8473742
## 52 35 - 39 0.8467756
## 53 35 - 39 0.9053645
## 54 35 - 39 0.7187136
## 55 35 - 39 0.7800328
## 56 35 - 39 0.8415641
## 57 40 - 44 0.9358703
## 58 40 - 44 0.8023624
## 59 40 - 44 0.8473742
## 60 40 - 44 0.8467756
## 61 40 - 44 0.9053645
## 62 40 - 44 0.7187136
## 63 40 - 44 0.7800328
## 64 40 - 44 0.8415641
## 65 45 - 49 0.9358703
## 66 45 - 49 0.8023624
## 67 45 - 49 0.8473742
## 68 45 - 49 0.8467756
## 69 45 - 49 0.9053645
## 70 45 - 49 0.7187136
## 71 45 - 49 0.7800328
## 72 45 - 49 0.8415641
## 73 5 - 9 0.9358703
## 74 5 - 9 0.8023624
## 75 5 - 9 0.8473742
## 76 5 - 9 0.8467756
## 77 5 - 9 0.9053645
## 78 5 - 9 0.7187136
## 79 5 - 9 0.7800328
## 80 5 - 9 0.8415641
## 81 50 - 54 0.9358703
## 82 50 - 54 0.8023624
## 83 50 - 54 0.8473742
## 84 50 - 54 0.8467756
## 85 50 - 54 0.9053645
## 86 50 - 54 0.7187136
## 87 50 - 54 0.7800328
## 88 50 - 54 0.8415641
## 89 55 - 59 0.9358703
## 90 55 - 59 0.8023624
## 91 55 - 59 0.8473742
## 92 55 - 59 0.8467756
## 93 55 - 59 0.9053645
## 94 55 - 59 0.7187136
## 95 55 - 59 0.7800328
## 96 55 - 59 0.8415641
## 97 60 - 64 0.9358703
## 98 60 - 64 0.8023624
## 99 60 - 64 0.8473742
## 100 60 - 64 0.8467756
## 101 60 - 64 0.9053645
## 102 60 - 64 0.7187136
## 103 60 - 64 0.7800328
## 104 60 - 64 0.8415641
## 105 65 - 69 0.9358703
## 106 65 - 69 0.8023624
## 107 65 - 69 0.8473742
## 108 65 - 69 0.8467756
## 109 65 - 69 0.9053645
## 110 65 - 69 0.7187136
## 111 65 - 69 0.7800328
## 112 65 - 69 0.8415641
## 113 70 - 74 0.9358703
## 114 70 - 74 0.8023624
## 115 70 - 74 0.8473742
## 116 70 - 74 0.8467756
## 117 70 - 74 0.9053645
## 118 70 - 74 0.7187136
## 119 70 - 74 0.7800328
## 120 70 - 74 0.8415641
## 121 75 - 79 0.9358703
## 122 75 - 79 0.8023624
## 123 75 - 79 0.8473742
## 124 75 - 79 0.8467756
## 125 75 - 79 0.9053645
## 126 75 - 79 0.7187136
## 127 75 - 79 0.7800328
## 128 75 - 79 0.8415641
## 129 80 - 84 0.9358703
## 130 80 - 84 0.8023624
## 131 80 - 84 0.8473742
## 132 80 - 84 0.8467756
## 133 80 - 84 0.9053645
## 134 80 - 84 0.7187136
## 135 80 - 84 0.7800328
## 136 80 - 84 0.8415641
## 137 85+ 0.9358703
## 138 85+ 0.8023624
## 139 85+ 0.8473742
## 140 85+ 0.8467756
## 141 85+ 0.9053645
## 142 85+ 0.7187136
## 143 85+ 0.7800328
## 144 85+ 0.8415641
## 145 Undefined 0.9358703
## 146 Undefined 0.8023624
## 147 Undefined 0.8473742
## 148 Undefined 0.8467756
## 149 Undefined 0.9053645
## 150 Undefined 0.7187136
## 151 Undefined 0.7800328
## 152 Undefined 0.8415641###################
# Graphs
trim <- function(ova) return(oapply(ova, NULL, mean, "Iter")@output)
ggplot(amount@output, aes(x=Gender, weight=amountResult, fill=Kala))+geom_bar()+
labs(
title="Kalansyönti Suomessa ikäryhmittäin",
y="Syönti (g/d)"
)
ggsave("Kalansyönti Suomessa ikäryhmittäin.svg")## Saving 7 x 5 in imageplot_ly(trim(total_amount), x=~Kala, y=~total_amountResult, color=~Kala, type="bar") %>%
layout(yaxis=list(title="Kalan kokonaiskulutus Suomessa (milj kg /a)"), barmode="stack")## Warning: `arrange_()` is deprecated as of dplyr 0.7.0.
## Please use `arrange()` instead.
## See vignette('programming') for more help
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_warnings()` to see where this warning was generated.plot_ly(trim(conc_vit), x=~Nutrient, y=~conc_vitResult, color=~Kala, type="scatter", mode="markers") %>%
layout(yaxis=list(title="Concentrations of nutrients (mg or ug /g)"))## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colorsggplot(conc_pfas@output, aes(x=Fish, y=conc_pfasResult))+geom_point()+ # Inputed data for missing species.
coord_flip()
ggplot(conc_pfas@output, aes(x=conc_pfasResult, color=Compound, linetype=Area))+stat_ecdf()+
scale_x_log10()+
stat_ecdf(data=conc_eukalat@output, aes(x=conc_eukalatResult))+
scale_linetype_manual(values=c("dotted","solid","twodash"))+
labs(
title="PFAS concentration in fishes in Finland",
x="PFAS concentration (ng/g fresh weight)",
y="Cumulative probability"
)
# The code may produce some negative values, which are removed from the graph
ggsave("PFAS-pitoisuus kalassa Suomessa.svg")## Saving 7 x 5 in imageggplot(conc@output, aes(x=concResult, colour=Fish))+stat_ecdf()+
facet_wrap(~Compound, scales="free_x")
#ggplot(oapply(expo_dir, cols=c("Iter"),FUN=mean)@output,
# aes(x=Gender, weight=expo_dirResult,fill=Fish))+geom_bar()+
# facet_wrap(~Compound, scales="free_y")+
# labs(title="Eri yhdisteiden saanti kalasta")
#ggsave("Yhdisteiden saanti kalasta Suomessa.svg")
plot_ly(trim(exposure), x=~Gender, y=~exposureResult, color=~Exposure_agent, text=~Exposure_agent, type="bar") %>%
layout(yaxis=list(title="Exposure to nutrients (g or ug /d)"))ggplot(trim(exposure), aes(x=Gender, weight=exposureResult, fill=Gender))+geom_bar()+
facet_wrap(~Exposure_agent, scales="free_y")+
labs(
title="Exposure to compounds",
y="(omega: mg/d; vit D: ug/d, PFAS: ng/d)"
)
cat("Kalaperäisiä tautitaakkoja Suomessa\n")## Kalaperäisiä tautitaakkoja Suomessaif(openv$N>1) {
tmp <- summary(oapply(BoDattr,NULL,sum,c("Age","Gender","Response")))
tmp <- data.frame(
Altiste = tmp$Exposure_agent,
Vaikutus = tmp$Response,
Keskiarvo = as.character(signif(tmp$mean,2)),
"95 luottamusväli" = paste0(signif(tmp$Q0.025,2)," - ", signif(tmp$Q0.975,2)),
Keskihajonta = signif(tmp$sd,2)
)#[rev(match(lev, tmp$Exposure_agent)),]
oprint(tmp)
tmp <- summary(oapply(BoDattr,NULL,sum,c("Age","Gender","Exposure_agent")))
tmp <- data.frame(
Terveysvaikutus = tmp$Response,
Keskiarvo = signif(tmp$mean,2),
"95 luottamusväli" = paste0(signif(tmp$Q0.025,2)," - ", signif(tmp$Q0.975,2)),
Keskihajonta = signif(tmp$sd,2)
)
oprint(tmp)
}
ggplot(trim(BoDattr), aes(x=Exposure_agent, weight=BoDattrResult, fill=Response))+geom_bar()
ggplot(trim(BoDattr), aes(x=Response, weight=BoDattrResult, fill=Exposure_agent))+geom_bar()
ggplot(trim(BoDattr), aes(x=Age, weight=BoDattrResult, fill=Response))+geom_bar(position="stack")
plot_ly(trim(BoDattr), x=~Exposure_agent, y=~BoDattrResult, color=~Response, text=~paste(Age, Exposure_agent, sep=": "), type="bar") %>%
layout(yaxis=list(title="Disease burden (DALY /a); CHD2=coronary heart disease"), barmode="stack")## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors################ Insight network
gr <- scrape(type="assessment")
objects.latest("Op_en3861", "makeGraph") # [[Insight network]]## Loading objects:
## makeGraphgr <- makeGraph(gr)## Loading required package: DiagrammeR## Loading objects:
## formatted
## Loading objects:
## chooseGr#export_graph(gr, "ruori.svg")
#render_graph(gr) # Does not work: Error in generate_dot(graph) : object 'attribute' not found
##################### Diagnostics
objects.latest("Op_en6007", code_name="diagnostics")## Loading objects:
## showind
## binoptest
## showLoctable
## ovashapetestshowLoctable()showind()## Hg is not an ovariable.
## subgrouping is not an ovariable.
## sumExposcen is not an ovariable.
## mc2d is not an ovariable.
## mc2dparam is not an ovariable.
## TEFversion is not an ovariable.